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Abstract 

In this article, we perform an extensive study of flavor observables in a two-Higgs-doublet model 

. (2HDM) with generic Yukawa structure (of type III). This model is interesting not only because it 

is the decoupling limit of the Minimal Supersymmetric Standard Model (MSSM) but also because 

of its rich flavor phenomenology which also allows for sizable effects not only in FCNC processes but 

also in tauonic B decays. We examine the possible effects in flavor physics and constrain the model 

both from tree-level processes and from loop-observables. The free parameters of the model are 

the heavy Higgs mass, tan/3 (the ratio of vacuum expectation values) and the " non-holomorphic" 

Qh Yukawa couplings (/ = u,d,£). In our analysis we constrain the elements in various ways: 

Q-r In a first step we give order of magnitude constraints on e{j from 't Hooft's naturalness criterion, 

finding that all efj must be rather small unless the third generation is involved. In a second 

step, we constrain the Yukawa structure of the type-Ill 2HDM from tree- level FCNC processes 

! (B s ,d ~^ A* - ') Kl — > -D — > fi + fJ-~, AF = 2 processes, r~ — > u~u + fi~, r~ — > e~fi + n~ 

and jjT — > e~e + e~) and observe that all flavor off-diagonal elements of these couplings, except 

e 32 31 an d e 23 13 mus t be very small in order to satisfy the current experimental bounds. In a 

third step, we consider Higgs mediated loop contributions to FCNC processes (b — > s(d)~f, B s ^ 

mixing, K — K mixing and u — > ej) finding that also e" 3 and e^ 3 must be very small, while the 

bounds on e 31 and €32 are especially weak. Furthermore, considering the constraints from electric 

dipole moments (EDMs) we obtain constrains on some parameters e"' . Taking into account the 

^ . constraints from FCNC processes we study the size of possible effects in the tauonic B decays 

5—i < 

(B — > tv , B — > Dtv and B — > D*tv ) as well as in — > tv, — > uv, K(ir) — > ev, K(ir) — > fiv 
and r —> K^)^ which are all sensitive to tree-level charged Higgs exchange. Interestingly, the 
unconstrained £32 31 are just the elements which directly enter the branching ratios for B —> tv , 
B — > Dtv and B — > D*tv . We show that they can explain the deviations from the SM predictions 
in these processes without fine tuning. Furthermore, B — >■ tv , B — > Dtv and B — > D*tv can even 
be explained simultaneously. Finally, we give upper limits on the branching ratios of the lepton 
flavor-violating neutral B meson decays (B Sj d — > fie, B s ^ — > re and B s c i — > tu) and correlate the 
radiative lepton decays (r — > /X7, r — > e^f and fi — > e^) to the corresponding neutral current lepton 
decays (r~ — > n~u + fi~, t~ — > e~fi + fi~ and pT — > e~e + e~). A detailed appendix contains all 
relevant information for the considered processes for general scalar-fermion-fermion couplings. 
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I. INTRODUCTION 



Two-Higgs-doublet models (2HDMs) [l| have been under intensive investigation for a 
long time (see for example Ref. {2] for an introduction or Ref. jsj for a recent review article). 
There are several reasons for this great interest in 2HDMs: Firstly, 2HDMs are very simple 
extensions of the Standard Model (SM) obtained by just adding an additional scalar SU{2)l 
doublet to the SM particle content. This limits the number of new degrees of freedom and 
makes the model rather predictive. Secondly, motivation for 2HDMs comes from axion 
models because a possible CP-violating term in the QCD Lagrangian can be rotated 
away Q if the Lagrangian has a global U(l) symmetry which is only possible if there are 
two Higgs doublets. Also the generation of the baryon asymmetry of the Universe motivates 
the introduction of a second Higgs doublet because in this way the amount of CP violation 
can be large enough to accommodate for this asymmetry, while the CP violation in the SM 
is too small ((J. Finally, probably the best motivation for studying 2HDMs is the Minimal 
Supersymmetric Standard Model (MSSM) where supersymmetry enforces the introduction 
of a second Higgs doublet [?| due to the holomorphic superpotential. Furthermore, the 
2HDM of type III is also the effective theory obtained by integrating out all super-partners 
of the SM-like particles (the SM fermion, the gauge boson and the Higgs particles of the 
2HDM) from MSSM. 

2HDMs are not only interesting for direct searches for additional Higgs bosons at collid- 
ers. In addition to these high energy searches at the LHC also low-energy precision flavor 
observables provide a complementary window to physics beyond the SM, i.e. to the 2HDMs. 
In this respect, FCNC processes, e.g. neutral meson decays to muon pairs (B s ^) — > /U'V - , 
D — )■ ji + [f and K L — > fi + fi~) are especially interesting because they are very sensitive to fla- 
vor changing neutral Higgs couplings. However, also charged current processes like tauonic 
5-meson decays are affected by the charged Higgs boson and b — > 57 provides currently the 
best lower limit on the charged Higgs mass in the 2HDM of type II. 

Recently, tauonic B decays received special attention because the BABAR collaboration 
performed an analysis of the semileptonic B decays B — > Dtv and B — > D*tu reporting a 
discrepancy of 2.0 a and 2.7 a from the SM expectation, respectively. The measurements of 
both decays exceed the SM predictions, and combining them gives a 3.4 a deviation from 
the SM [8|, |9| expectation, which constitutes first evidence for new physics in semileptonic 
B decays to tau leptons. This evidence for the violation of lepton flavor universality is 
further supported by the measurement of B — > rv by BABAR 10, 11 ] and BELLE [12I, LX3| 



which exceeds the SM prediction by 1.6 a using V uh from the global fit [14J]. Assuming that 
these deviations from the SM are not statistical fluctuations or underestimated theoretical or 
systematic uncertainties, it is interesting to ask which model of new physics can explain the 
measured values. Since, a 2HDM of type II cannot explain B — > tv , B — > Dtv and B — > 
D*tu simultaneously [8[, one must look at 2HDMs with more general Yukawa structures. 
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Also 2HDMs of type III with Minimal Flavor Violation (MFV) 15[ cannot explain these 
deviations from the SM but a 2HDM of type III (where both Higgs doublets couple to 
up quarks and down quarks as well) with flavor-violation in the up sector, is capable of 
explaining B — > tv , B — > Drv and B — > D*tv without fine tuning 16]. 

These points motivate us to perform a complete analysis of flavor-violation in 2HDMs 
of type III in this article. For this purpose we take into account all relevant constraints 
from FCNC processes (both from tree-level contributions and from loop-induced effects) 
and consider afterwards the possible effects in charged current processes. 

This article is structured as follows: In Sec. [Til we review the Yukawa Lagrangian of the 
2HDM of type III. In Sec. IHII we give a general overview on the constraints on 2HDMs and 
update the bounds on the 2HDM of type II. The following sections discuss in detail the 
constraints on the 2HDM of type III parameter space from 't Hooft's naturalness argument 
(Sec. II Vj) . from tree-level FCNC processes (Sec. |V|) and from loop-induced charged and 
neutral Higgs mediated contributions to the flavor observables (Sec. IVIj) . Sec. IVHI studies 
the possible effects in charged current decays [B — > tv , B — > Dtv , B — \ D*tv , Di s \ — > tv, 
D( s ) fiv, K(ti) — > ev, K{j) — > fiv, r — > K{ji)v) and Sec. IVIIII is denoted to the study 
of the upper limits on the branching ratios B St d — > Tfi, B St d — > re, B St d — >■ fie and the 
correlations among r~ — > — > e~fi + fi~, fi~ — > e~e + e~ and r — > fij, t — > ej, 
fi — > e'y. Finally, we conclude. A detailed appendix contains some of the input parameters 
used in our analysis, general expressions for some branching ratios as well as all the relevant 
Wilson coefficients for b — > s(d)'j, AF = 2 processes, leptonic neutral meson decays (AF = 
1), LFV transitions, EDMs, anamolous magnetic moment (AMM) of muon and (semi-) 
leptonic charged meson decays for general charged and/or neutral scalar- fermion-fermion 
couplings. 



II. SETUP 

The SM contains only one scalar weak-isospin doublet, the Higgs doublet. After elec- 
troweak symmetry breaking its vacuum expectation value ("vev") gives masses to up quarks, 
down quarks and charged leptons. The charged (CP-odd neutral) component of this dou- 
blet becomes the longitudinal component of the W (Z) boson, and thus we have only one 
physical CP-even neutral Higgs particle in the SM. In a 2HDM we introduce a second Higgs 
doublet and obtain four additional physical Higgs particles (in the case of a CP conserving 
Higgs potential): the neutral heavy CP-even Higgs H°, a neutral CP-odd Higgs A and the 
two charged Higgses H ± . 

As outlined in the introduction we consider a 2HDM with generic Yukawa structure 
(2HDM of type III). One motivation is that a 2HDM with natural flavor-conservation (like 
type I or type II) cannot explain B — > Dtv , B — > D*tv and B ^ tv simultaneously, while 



the type III model is capable of doing this [161 ] . Beside this, our calculations in the 2HDM 
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Ill are the most general ones in the sense that they can be applied to models with specific 



flavor-structures like 2HDMs with MFV 



15. 



3,13]. 



In this sense also our bounds are model 



independent, because they apply to any 2HDM with specific Yukawa structures as well (in 
the absence of large cancellations which are unlikely). Finally the type-Ill 2HDM is the 
decoupling limit of the MSSM and the calculated bounds can be translated to limits on the 
MSSM parameter space. 

The fact that the 2HDM III is the decoupling limit of the MSSM also motivates us to 
choose for definiteness a MSSM like Higgs potential 1 which automatically avoids dangerous 



CP violation. The matching of the MSSM on the 2HDM Yukawa _sector has 
in detail. For the MSSM with MFV it was calculated in Ref. 



19 



Deen considered 



21 



22 



23 



for the MSSM with generic flavor structure in Ref. 25[ (neglecting the effects of the A 



24J and 



terms) and in Ref. 26[ (including the A-terms). 
were calculated for the flavor- conserving case in 



ven the next-to-leading order corrections 



general MSSM inRef. 
considered 1291 



30 



found to be small 



31 



3s 



28 
32 



271 ] and for the flavor- changing one in the 



Also the one-loop corrections to the Higgs potential have been 



33 



34 



35 



36 



37J, but their effects on flavor-observables were 



Following the notation of Ref. 26|, |28|, [39[ we have the following Yukawa Lagrangian in 



the 2HDM of type III starting in an electroweak basis: 

c Y = Q a fL [Yfr^Hf - ejrK] d iR + Q°j L [Yfr^ ab H b : - e u f rm\ UiR + h. c . . (i) 

Here a, b denote SU{2) L - indices, e a & is the two-dimensional antisymmetric tensor with 
ei2 = — 1 and the Higgs doublets are defined as : 





-\ 


'hi 










\- 


(* 






{ H° 

\ u 



with (Hd) 



with (H u ) 




(2) 



Apart from the holomorphic Yukawa-couplings Y^ ew and V^ ew , we included the non- 
holomorphic couplings e^ ew (q = u, d) as well. 

As a next step we decompose the SU (2) doublets into their components and switch to a 



1 If wc would require that the Higgs potential possesses a Z2 symmetry the results would be very similar. 
The heavy Higgs masses squared would still differ by terms of the order of v 2 and only Higgs self-couplings 
would be different, but they do not enter the flavor-processes at the loop-level under consideration. 
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basis in which the holomorphic Yukawa couplings are diagonal: 



C Y = -d JL [Y*6 fi H* + ~e% Hi] d lR - u fL [Y^6 fi H°* + eft H°] u iR 
+ Z fL V, j [Y*5 ji -cotl3Zj i ]H?d iR 
+ d fL V* f [Y^Sji - tan/3e£] H^u iR + h.c. . 



(3) 



where, tan/3 = v u /vd is the ratio of the vacuum expectation values v u and Vd acquired 
by H u and Ha, respectively. We perform this intermediate step, because this is the basis 
which corresponds to the super-CKM basis of the MSSM and the couplings efj can be directly 
related to loop-induced non-holomorphic Higgs coupling. The wave-function rotations U \ ' 
necessary to arrive at the physical basis with diagonal quark mass matrices are defined by 



U jf L * m jk U ki R =™jfi- 
They modify the Yukawa Lagrangian as follows: 



(4) 



Cy — — df l 

+ u fL V fj 

+ d fL v; f 



5 fl -e a fl t^P)Hf + e%H° u 



di 



Vd 

m 



Sji - (cot (3 + tan 0) 



v 



-S^- (tan/3 + cot/3)e" 



Hd* d iR 

Hl*u iR + h.c. . 



Here, m qi are the physical running quark masses and 



(5) 



V fi = VffUjf , (6) 

is the CKM matrix. The Higgs doublets H u and project onto the physical mass eigen- 
states H° (heavy CP-even Higgs), h° (light CP-even Higgs), A (CP-odd Higgs) and H ± in 
the following way: 

if ° = —= (H° sin a + h° cos a + iA° cos (3) , 
v2 

H° d = — (H° cos a - h° sin a + iA° sin (3) , 
V2 



Hi = cos/3if H 
Hj = sin /3 H~ 



(7) 



where, a is the mixing angle necessary to diagonalize the neutral CP-even Higgs mass matrix 
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(see e.g. 



40|). Since we assume a MSSM-like Higgs potential 2 we have 



4- a Vu 
tan p — — , 

tan 2a = tan 2/3 — ^ f , ( 8 ) 

m iv± = 171 a° + M w » = m^o + M| - m 2 , 

, — 7T „ 7T 

with < a < and < 8 < —. 

2 y 2 

This means that in the phenomenologically interesting and viable limit of large val- 
ues of tan/3 and v <C m^o we have to a good approximation 3 : 

tan/3 ps — cot a , ^ 
m#o ps m#± ps m^o = m# . 

Without the non-holomorphic corrections efy, the rotation matrices Uj^' R would simulta- 
neously diagonalize the mass terms and the neutral Higgs couplings in Eq. (JSJ). However, in 
the presence of non-holomorphic corrections, this is no longer the case and flavor changing 
neutral Higgs couplings are present in the basis in which the physical quark mass matrices 
are diagonal. 

The Yukawa Lagrangian in Eq. fl5]) leads to the following Feynman rules 4 for Higgs-quark- 
quark couplings 



Hr™»P R + r%»p L ) (io) 



with 



J=l v Vd J 

T L d f u f = £cos/3 Vft - e>n/?) . (11) 

7=1 \ u 



.7 

Similarly, for the lepton case, the non-vanishing effective Higgs vertices are 



2 MSSM-like Higgs potential implies that in the large tan/3 limit and for v <C mn the charged Higgs mass 
mjj± , the heavy CP even Higgs mass m H a and the CP odd Higgs mass m^o are equal. 

3 For the SM-like Higgs boson h° we use m^o w 125 GeV in our numerical analysis. 

4 Hcrmiticity of the Lagrangian implies the relation Fff^* 1 = ^q^f 1 *- 



s 



(12) 



= E sin/3 U™ s ( ^ - 4, tan 

jr = l V V d 

Here, = (if , h°, A ) and the coefficients x k q are given by 

fc ( 1 • 1 * a\ 

x,. = =sina, 7=cosa, —=cosp , 

V y/2 y/2 V2 J n ,x 

* ( 1 1 • * • R 

Xi = 7=cosa, — =sma, — =sinp 

V y/2 V2 y/2 

This means that flavor-violation (beyond the one already present in the 2HDM of type II) 
is entirely governed by the couplings e\j . If one wants to make the connection to the MSSM, 
the parameters e\j will depend only on SUSY breaking parameters and tan/3. 

III. CONSTRAINTS ON THE 2HDM PARAMETER SPACE 
- GENERAL DISCUSSION AND OVERVIEW 

In this section we give an overview on flavor observables sensitive to charged Higgs con- 
tributions. We review the constraints on the 2HDM of type II and discuss to which extent 
these bounds will hold in the 2HDM of type III. A detailed analysis of flavor constraints on 
the type-Ill 2HDM parameter space will be given in the following sections. 

The most common version of 2HDMs, concerning its Yukawa sector, is the 2HDM of type 



II which respects natural flavor conservation 4l|] by requiring that one Higgs doublet couples 



only to up-quarks while the other one gives masses to down-type quarks and charged leptons 
(like the MSSM at tree- level). Flavor-observables in 2HDMs of type II have been studied 
in detail 42j, |43|, |44(. In the type II model there are no tree-level flavor- changing neutral 
currents and all flavor violation is induced by the CKM matrix entering the charged Higgs 
vertex. In this way the constraints from FCNC processes can be partially avoided. This is 
true for AF = 2 processes where the charged Higgs contribution is small, for K L — > 
D° — > fi + fi~ (due to the tiny Higgs couplings to light quarks) and all flavor observables in 
the lepton sector. However, the FCNC processes b — > 57 (also to less extent b — > dj) and 
B s — > / u + fi~ are sensitive the charged Higgs contributions. In addition, direct searches at 
the LHC and charged current processes restrict the type-II 2HDM parameter space. 

Among the FCNC processes, the constraints from 6—7-57 are most stringent due to the 



necessarily constructive interference with the SM contribution [45l , l46l . 1471 . |48I |. The most 
recent lower bound on the charged Higgs mass obtained in Ref. 49J is m#± > 360 GeV 
which includes NNLO QCD corrections and is rather independent of tan/3. In the type- 
Ill 2HDM this lower bound on the charged Higgs mass can be weakened due to destructive 
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interference with contributions involving e?-. Also in B s — > fi + fi~ (and B^ — >■ a sizable 

loop-induced effect is possible in the 2HDM II, but the constrains are still not very stringent 
even if the new LHCb measurement are used. The reason for this is that, taking into account 
the constraints from b — > S7 on the charged Higgs mass, the branching ratio for B s — »■ 
in the 2HDM II is even below the SM expectation for larger values of tan (3 50|, |5ll [52J due 
to the destructive interference between the charged Higgs and the SM contribution. 

Regarding charged current processes, tauonic B decays are currently most sensitive to 
charged Higgs effects. Here, the charged-Higgs contribution in the type-II 2HDM to B — > 



rv interferes destructively with the SM contribution 53|, |54|. The same is true for B 



D*TV 



551 ] and B Dtv 



42 



56 



57|. As outlined in the introduction this leads to the 



act that the 2HDM II cannot explain B — > tv , B — > Dtv and B — > D*tv simultaneously 
8|. Other charged current observables sensitive to charged Higgses are D( s ) — > fiv, — > 



TV 



58 



601 ] . r —> K{n)v and K — > jiv/n — > \iv [6l| (see (44j for a global analysis) 



2HDM II 



100 



80 



_ 60 
^ 40 



20 





200 400 600 800 1000 
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FIG. 1: Updated constraints on the 2HDM of type II parameter space. The regions compatible 
with experiment are shown: b — > S'j (yellow), B — > Dtv (green), B — )• tv (red), B s — > fi + ' [i~ 
(orange), K — >■ [iv / ir — )• fiv (blue) and B — > D*tv (black). Note that no region in parameter space 
is compatible with all processes. Explaining B — > D*tv would require very small Higgs masses 
and large values of tan (3 which is not compatible with the other observables. To obtain this plot, 
we added the theoretical uncertainty linear on the top of the 2 a experimental error. 




Fig. [TJ shows our updated constraints on the 2HDM II parameters space from b — > sj, 
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B — > rv , B — > Dtv , B —¥ D*tv , B s — > fi + fi~ and K — > fiv/n — > \iv. We see that in 
order to get agreement within 2 a between the theory prediction and the measurement of 
B — > D*tv , large values of tan (3 and light Higgs masses would be required which is in 
conflict with all other processes under consideration. 

Concerning direct searches the bounds on the charged Higgs mass are rather weak due to 
the large background from W events. The search for neutral Higgs bosons is easier and the 
CMS bounds 5 on m^o from A — > t + t~ are shown in Fig. [2j These bounds were obtained 
in the MSSM, but since the MSSM corrections to A — > t + t~ are rather small and since we 
consider a MSSM-like Higgs potential, these bounds also hold in the 2HDM III as long as 
the Peccei-Quinn symmetry breaking in the lepton sector is small. 




FIG. 2: Plot from the CMS collaboration taken from Ref. 62|: allowed regions in the m^o-tan/3 
plane from ^4° — > t + t~ . The analysis was done in the MSSM, but since we consider a 2HDM with 
MSSM-like Higgs potential and the MSSM corrections to the A°tt vertex are small, we can apply 
this bound to our model. However, a large value of e| 3 in the 2HDM of type III could affect the 
conclusions. 



Going beyond the simp 



e Yukawa structure of the 2HDM of type II, also 2HDMs of 
type III with MFV 15|, |l7|, ll8|, alignment 63|, |64( or natural flavor conservation 
have been analyzed in detail. However, flavor-observables in type III models with generic 



Note that we did not use the bounds from unpublished CMS update of the ^4° — > t + t analysis. 
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flavor-structure have received much less attention. Ref. 65[ considered the possible effects 
of the flavor-diagonal terms and Ref. (ssj) considers leptonic observables. As outlined in 
the introduction, 2HDMs of type II (or type III with MFV) cannot explain B — > Dru and 
B — > D*tu simultaneously jsj (and for B — > tv fine tuning is needed fla]). 

In the following sections we will study in detail the flavor-observables in the 2HDM with 
generic flavor-structure 67[, but for definiteness, with MSSM-like Higgs potential. For this 
purpose, all processes described above are relevant. In addition, AF = 2 processes, lepton 
flavor violating observables (LFV), EDMs, r — > K(tt)u/ K(tt) — > fiu and K — > [i[e)vj-n — > 
\i{e)v will turn out to give information on the flavor structure of the 2HDM of type III. 
Furthermore, we will investigate to which extent contributions to B Sidi — > Tfi, B Sfd — > re, 
B St d A*e and muon anomalous magnetic moment are possible. 



IV. CONSTRAINTS FROM 'T HOOFT'S NATURALNESS CRITERION 

The naturalness criterion of 't Hooft states that the smallness of a quantity is only 
natural if a symmetry is gained in the limit in which this quantity is zero. This means on 
the other hand that large accidental cancellations, which are not enforced by a symmetry, are 
unnatural and thus not desirable. Let us apply this reasoning to the fermion mass matrices 
in the 2HDM. We recall from the last section the expressions for the fermion mass matrices 
in the electroweak basis: 



m % 


= v d Yf™ 


1 dew 




— „, -yuew 

— Vulij 


+ v d e i:j 




= v d Y^ 


+ V f iew 



Diagonalizing these fermion mass matrices gives the physical fermion masses and the 
CKM matrix. Using 't Hooft's naturalness criterion we can demand the absence of fine- 
tuned cancellations between v d Y^' {v u Y^) and v u efj (vaefj). Thus, we require that the 
contributions of v u eff and i^e"- to the fermion masses and CKM matrix do not exceed the 
physical measured quantities. 

In first order of a perturbative diagonalization of the fermion mass matrices, the diagonal 
elements give rise to the fermion masses, while (in our conventions) the elements m{- 
with i < j (i > j) affect the left-handed (right-handed) rotations necessary to diagonalize 
the fermion mass matrices. The left-handed rotations of the quark fields are linked to the 
CKM matrix and can therefore be constrained by demanding that the physical CKM matrix 
is generated without a significant degree of fine-tuning. However, the right-handed rotations 
of the quarks are not known and the mixing angles of the PMNS matrix are big so that 
for these two cases we can only demand that the fermion masses are generated without too 
large accidental cancellations. Note, that in Eq. f)14p the elements e{ J ew enter, while the 
elements e{ which we want to constrain from flavor observables are given in the physical 
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basis with diagonal fermion masses. This means that in order to constrain from 't Hooft's 
naturalness criterion we have to assume in addition that no accidental cancellation occur by 
switching between the electroweak basis and the physical basis. In conclusion this leads to 
the following upper bounds 



d(u) 



< |Vj KM | x max [m di{Ui) ,m d]{U:j) ] for % < j 

< max [m di{Ui)l m d]{Uj) ] for i>j, 
u4j \ < max [m^m^] . 



d(u) 
V u(d) e ij 



(15) 



In the large tan/3 limit, inserting the quark masses m q (fi) at the Higgs scale (which we 



choose here to be ^m g g 
from Eq. flTSJ): 



500 GeV), we can immediately read off the upper bounds on ' 



/l.3xl0~ 4 5. 



- 1 j I 



< 



x 10" 



5.1 x 10^ 5 \ 



2.6 x 10~ 4 2.6 x 10~ 4 5.9 x 10~ 4 
\l.4 x 10~ 2 1.4 xlO^ 2 1.4x10-2/.. 



/3.4xl0~ 4 3.2 x 10~ 2 1.6xl0 _1 \ 



- 'j 



< (tan/3/50) 



1.4 x 10" 1 1.4 x 10" 1 



1.9 



V 



(16) 



< 



/2.9xl0~ 6 6.1 x 10^ 4 1.0 x 10~ 2 \ 

6.1 x 10~ 4 6.1 x 10^ 4 1.0 x 10~ 2 
\1.0 x 10~ 2 1.0 x 10~ 2 1.0 x 10~ 2 y 



Of course, these constraints are not strict bounds in the sense that they must be respected 
in any viable model. Anyway, big violation of naturalness is not desirable and Eq. (fT6|) gives 
us a first glance on the possible structure of the elements e{- 
possible to explain B — > rv , B — > Drv and B — > D*tv using e 
Eq. (TT6]) . while if one wants to explain B — >■ rv with 't Hooft's naturalness criterion is 
violated. 



As we will see later, it is 
3i 32 without violating 



V. CONSTRAINTS FROM TREE-LEVEL NEUTRAL-CURRENT PROCESSES 

The flavor off-diagonal elements e|j- (with i / j) give rise to flavor- changing neutral 
currents (FCNCs) already at the tree-level. Comparing the Higgs contributions to the 
loop-suppressed SM contributions, large effects are in principle possible. However, all ex- 
perimental results are in very good agreement with SM predictions, which put extremely 
stringent constraints on the non-holomorphic terms ef-. 
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In this section we consider three different kinds of processes: 

• Muonic decays of neutral mesons (B sd — > fi + fi~, K^ —> and D° — > fi + fi~). 

• AF = 2 processes (D — D, K — K , B s — B s and B d — B d mixing). 

• Flavor changing lepton decays (r~ — )■ r~ — > e~fi + fi~ and yT — > e~e + e~). 

As we will see in detail in Sec. IV A\ the leptonic neutral meson decays B s d — > 

and D° — >■ /i + fi put constraints on the elements efj (with i j) and £"2^1 
already if one of these elements is non-zero, while B d —B d , B s —B s , K—K and D—D mixing 
only provide constraints on the products e^e^* and £i 2 e 2i (Sec. IV Bp . This means that the 
constraints on AF = 2 processes can be avoided if one element of the product e^ej* is zero, 
while the constraints from the leptonic neutral meson decays can only be avoided if the 
Peccei Quinn symmetry breaking for the leptons is large such that ~ Tn^/vu is possible. 

In Sec. IV CI we will consider the flavor changing lepton decays t~ — > fi~fi + fi~, t~ — > 
e~fi + fi~ and fi" — > e~e + e~ which constrain the off-diagonal elements £2332; e i3 3i an d e i2 2i> 
respectively. 



A. Leptonic neutral meson decays: B s d — > /U + /i , Kl — > y + fi and D° — > fi + fi 



Muonic decays of neutral mesons {B s — > B d — > fi + fi~, K^ — > fi + fi~ and D° — > 

are strongly suppressed in the SM for three reasons: they are loop- induced, helicity 
suppressed and they involve small CKM elements. Therefore, their branching ratios (in 



the SM) are very small and in fact only Kl — > fi + n~ and recently also B s — > 68] 
have been measured, while for the other decays only upper limits on the branching ratios 
exist (see Tabled]). We do not consider decays to electrons (which are even stronger helicity 
suppressed) nor B ds — > t + t~ (where the tau leptons are difficult to reconstruct) because the 
experimental limits are even weaker. The study of meson decays to lepton flavor-violating 
final states is postponed to Sec. IVIIH 

We see from Fig. [3] that the off-diagonal elements of ef 3 31 , £33 32 , £12 21 anc ^ e i2 21 directly 
give rise to tree- level neutral Higgs contributions to B d — > B s — > , K L — >• 

and D° — > fi + fi~, respectively. 

In principle, the constraints from these processes could be weakened, or even avoided, if 
~ rne 2 / v u- Anyway, in this section we will assume that the Peccei Quinn breaking for 
the leptons is small and neglect the effect of q 2 in our numerical analysis for setting limits 



on e\ r 
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Process 


Experimental value 


SM prediction 


B [B s -> n+/i-] 


3.2±^ x 10- 9 [68] 


(3.23 ± 0.27) x 10- 9 [69J 


B [B d -> 


< 9.4 x 10~ 10 (95% CL) [68] 


(1.07 ±0.10) x 10~ 10 [69] 


B -> /i+Z^kort 


< 2.5 x 10" 9 [7pJ 


« 0.9 x 10" 9 [7pJ 


B [-D -»• 


< 1.4 x 10~ 7 (90% CL) [71] 





TABLE I: Experimental values and SM predictions for the branching ratios of neutral meson decays 
to muon pairs. For Kl — > [i + [i~ we only give the upper limit on the computable short distance 
contribution [70(| extracted from the experimental value (6.84 ± 0.11) x 10~ 9 (90% CL) [Tlf| . The 
SM prediction for D° — > ji + fj,~ cannot be reliably calculated due to hadronic uncertainties. 




FIG. 3: Feynman diagram showing the neutral Higgs contribution to B s ^ — > /i + ji , Kl — > 
and D° — > 



1. B s4 ->■ n+n 

For definiteness, consider the decay of a neutral B s (5s) meson (the corresponding decay 
of a Bd meson follow trivially by replacing s with d and 2 with 1) to a moun pair. The 
effective Hamiltonian governing this transition is 6 

f 2 M 2 

ft = L_W_ [cbsQbs + QbsQbs + C bsQbs + C'bsQ'bs + C'bsQlbs + Q'bs Q'^] + h.C. , 

(17) 



6 The complete expression for the Hamiltonian and the branching ratio including lepton flavor-violating 
final states is given in the appendix. 
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where the operators are defined as 



01 



(bP L s) {Jin) , 
(bP L s) (fiysfi) , 



(18) 



and the primed operators are obtained replacing Pl with Pp. The corresponding expression 
for the branching ratio in terms of the Wilson coefficients reads 



B [B s -»• n+n- 



G%Mf, 



4— 



57T U 



Ml 



h fB a m l T B S 



X 



Ml (C b P s - Cp) 



2 (m b + m s ) m h 



(jrbs^ 



MUC 



bs 



2 (m b + m s ) m M 



x 



2 

mi, 
4^ 



m 



B s 



(19) 

Concerning the running of the Wilson coefficients due to the strong interaction, the operators 
O b / and 0% s correspond to conserved vector currents with vanishing anomalous dimensions. 
This means that their Wilson coefficients are scale independent. The scalar and pseudo- 
scalar Wilson coefficients Cg s and Cp {C^ s and Cp s ) have the same anomalous dimension 
as quark masses in the SM which means that their scale dependence is given by: 



m q (ni 



high ) 



high ) 



(20) 



where m q is the running quark mass with the appropriate number of active flavors. In the 
SM, Ca is the only non-vanishing Wilson coefficient 



ribs 



V&VuY 



m t 



v; b v cs Y 



Ml 



(21) 



where, the function Y is defined as Y = tjyYq such that the NLO QCD effects are included 
in 7] Y = 1.0113 



and the one loop Inami-Lim function Yq reads 72] 



Y n (x) = - 



A - x 



3x 



1 — x ( 1 — x) 



ln(x) 



(22) 



The complete Wilson coefficients for general quark- quark- scalar couplings are given in the 
appendix. In the 2HDM of type III, in the case of large tan/3 and v <C m#, the terms 
involving e\- simplify to 



ribs 



c b P s 



7T 



1 m e3 



rubs rubs 



GlM^2m 2 H 



7T 



V 



t 23 



tan 2 /3 



(23) 



G F M^2m 2 H 



-^e^tan 2 /?. 
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B s ^>n + ji (tan/?=30) 



B s ^n + ji (tan/?=30) 



nr. 




Re[e3 2 ] x 1 5 
B s ^n + iT (tan/3=50) 




Re[ef 2 ] x 10 5 




Re[E2 3 ] x 1 5 
B s ^n + iT (tan/3=50) 




Re[e^ 3 ] x 10 5 



FIG. 4: Allowed regions in the complex e 23 32 -plane from B s — > n + n for tan /3 = 30, tan /3 = 50 
and m H = 700 GeV (yellow), m H = 500 GeV (red) and m H = 300 GeV (blue). Note that the 
allowed regions for e 32 -plane are not full circles because in this case a suppression of B [B s — > [i + fi~] 
below the experimental lower bound is possible. 



To these Wilson coefficients the well known loop-induced type II 2HDM contributions 7 



ribs 



C b P s 



m b V£Vts Try 
2 2M^ 



tan 2 (3 



log (m 2 H / m 



m 2 H /m 



1 



(24) 



have to be added as well 52[. Note that since we give the Wilson coefficients at the matching 
scale, also mj and m t must be evaluated at this scale. 

We can now constrain the elements e 23 32 and ef 3 31 by demanding that the experimental 
bounds are satisfied within two standard deviations for B s — > or equivalently at the 

95% CL concerning B d — > The results for the constraints on e 23 and ef 2 (ef 3 and e 31 ) 

from B s — > fi + fi~ (Bd — > /i + /i~) are shown in Fig. H] (Fig. [5]). 

All constraints on 6i 331 and e 23 32 are very stringent; of the order of 10~ 5 . Both an 
enhancement or a suppression of B [B ds — > fi + fi~] compared to the SM prediction is possible. 



7 Since we want to put constraints on the elements ef 3 2 z we assume that the loop-induced 2HDM II 
contribution is not changed by elements e" 3 or £33. 
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-2-10123 -2-10 1 2 

Re^JxIO 5 Re[ef 3 ]x10 5 



FIG. 5: Allowed regions in the complex ef 3 31 -plane from — > /i + /i for tan/3 = 30, tan/3 = 50 
and m H = 700 GeV(yellow), m H = 500 GeV(red) and m H = 300 GeV(blue). 

While in the 2HDM II the minimal value for B [B^ s is half the SM prediction, in 

the 2HDM III also a bigger suppression of Bd yS — > is possible if ef 3 23 ^ 0. In principle, 

the constraints on e 23 (ef 3 ) from B s ^ — > are not independent of e 32 (e 31 ). Anyway, 

in the next section it will turn out that the constraints from AF = 2 processes are more 
stringent if both e 32 and e 23 are different from zero (the same conclusions hold for e 31 13 , 
e 2i,i2 an d e 21jl2 ). 

B s — > and Bd — > fi + fi~ can also be used to constrain the leptonic parameter e 22 . 

We will discuss the corresponding subject in SecJVIl 



2. K L fjL+fi- 

Concerning — > fi + fi~, the branching ratio and the Wilson coefficients can be obtained 
by a simple replacement of indices from Eq. ( fl9|) . Eq. ( 12T1) and Eq. ( )23|) . Due to the presence 
of large non-perturbative QCD effects, we require that the 2HDM III contribution together 
with the short distance piece of the SM contribution does not exceed the upper limit on 



the short distance contribution to the branching ratio calculated in Ref. 70]. The resulting 



constraints on ef 2 21 are shown in Fig. [6j They are found to be extremely stringent (of the 
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-3-2-1012 -3-2-1012 
Re[ef 2 ]x10 6 Rele^lxlO 6 



FIG. 6: Allowed regions in the complex ef 2 2 i-pl an e from Kl — > fi + fi for tan/3 = 30, tan/3 = 50 
and m H = 700 GeV(yellow), m H = 500 GeV(red) and m H = 300 GeV(blue). 

order of 10 -6 ). 



3. D° fi + fjL~ 

The analogous expressions for the branching ratio for D° — > (D°(cu)) follow by a 

straightforward replacement of indices in Eq. ffT9]) but the Wilson coefficients in the type-Ill 
2HDM for D° — > have a different dependence on tan f3: 



2 l 
(j F 1V1 W Zm H (25) 



2 i t 

°s "°p - n 2 t\/t2 7T~2 ~ e 21 tanp. 



G 2 F Mw 2m 2 H 

Differently than for B^ s At 4 'fj,~ the SM contribution cannot be calculated due to non- 
perturbative effects and the 2HDM II contribution is numerically irrelevant. Since we do 
not know the SM contribution, we require that the 2HDM III contribution alone does not 
generate more than the experimental upper limit on this branching ratio. 
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r°,h°,A° 



- q * (f q \ 
if \ e fi) 



FIG. 8: Feynman diagram contributing to B^-B^^, K-K and D-D mixing. 



It is then easy to express the constraints on e\ 2 21 in terms of the parameters m# and 
tan (3: 

The resulting bounds on e" 2 21 (setting one of these elements to zero) are shown in Fig. [71 



B. Tree-level contributions to AF = 2 processes 

In the presence of non-zero elements e^- neutral Higgs mediated contributions to neutral 
meson mixing {Bd, s ~Bd,si K-K and D D mixing) arise (see Fig. [8]). In these processes, the 
2HDM contribution vanishes if the £7(1)pq symmetry is conserved. This has the consequence 
that the leading tan /3-enhanced tree-level contribution to the AF = 2 processes (shown 
in Fig. [8]) is only non-vanishing if e\- and are simultaneously different from zero (in 
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the approximation m^o = m#o and cot/3 = 0). Making use of the effective Hamiltonian 
defined in Eq. f[M|) of the appendix we get the following contributions to B s -B s mixing (the 
expressions for B^-Bd and A- A mixing again follow by a simple replacement of indices): 



C 4 



e 23 e 32 ,2 



m 



tan f3 . 



(27) 



H 



All other Wilson coefficients are sub-leading in tan (5. For D mixing, again only C 4 is non- 
zero and given by 

Cd = - e J^L . (28) 



m 



H 



After performing the renormalization group evolution |73|, [75! . l76l . l77j (here we used 
/jLh = 500 GeV at the high scale) it turns out that the dominant contribution to the hadronic 
matrix elements stems from 4 . Inserting the bag factors Tsl 
lattice QCD (see Table. IX Ij) . we get for the 2HDM of type III contribution 

(B° d \ C 4 A \Bf) w 0.26 C 4 GeV 3 , 
(B° s \ C A A \B° S ) w 0.37 C 4 GeV 3 , 
(K°\ C 4 4 I A' ) w 0.30 C 4 GeV 3 , 
(D°\C 4 4 \D°) ^0.18C 4 GeV 3 , 



79l | and decay constants from 



(29) 



where, we used the normalization of the meson states as defined for example in 761 ] . In 
Eq. (1291) the Wilson coefficients within the matrix elements are at the corresponding meson 
scale while C 4 on the right-handed side is given at the matching scale ma- For computing 
the constraints on e f 3 eo *, e^e^ and ef 2 e 2i we use ^he online update of the analysis of the 



UTfit collaboration 



C Bo e 



for Bd — Bd and B s — B s mixing and 



801] 8 . For this purpose we define 



2i<pi 



1 + 



[B° q \n!f f \B0) 



(30) 



a 



1 + 



c 



AM* 



Im [(K°\H? f p f \K° 
Im [{K \n s ef M f \K )] 
Re \(K°\ '-. 

1 + 



Re [(K° 


njNP 
H eff 




Re [(A0|^ 


A0>] 



(31) 



See also the online update of the CKMfitter group for an analogous analysis [14 1. 
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for K — K mixing. Using for the matrix elements of the SM Hamiltonian [81] 

(B° d \H^= 2 \B° d ) » (1.08+ 1.25z) x 10" 13 GeV, 

(B° s \n^r 2 \B° s ) « (59-2.2^) x 10~ 13 GeV, (32) 

(*. |?^= 2 |if ) ~ (115 + 1.16^ x 10~ 17 GeV, 

we can directly read off the bounds on C4 and thus on ef 2 e 2i> e i3 e 3i anc ^ e 23 e 32 : 

- 2.0 x 10- 10 < Re [441 ( tan/3/ f° ^ 2 < 6.0 x lO" 10 , (33) 

L 23 32j ^ mH / 500 GeV J - ' v ; 

-3.0X10- < [m [44] (_^L) 2 < 7.0 x 10- , (34) 

-3.0xl0-< Re [44] (^^) 2 < 1.5 x 1(T« , (35) 

-1.5xl0-< I m [44](-^^y<2. 5 xl0-", (36) 

-1.0X10- < R e[44](^^^) 2 <3.0xl0-, (37) 

-4.0xl0-< I m [44](^|^) 2 <2.5xl0-. (38) 

We see that if el is of the same order as el these bound are even more stringent than the 
ones from B ds — > and Kl — > fi + n~ computed in the last subsection. 

For D — D mixing, the SM predictions is not known due to very large hadronic uncer- 
tainties. In order to constrain the NP effects we demand the absence of fine tuning, which 
means that the NP contribution, which are calculable short distance contributions, should 
not exceed the measured values. Concerning the 2HDM III contribution, there is no tan /3 
enhancement and taking into account the recent analysis of UTfit collaboration 82| we arrive 
at the following constraints (for m# = 500 GeV): 



I e i2 e 2il < 2.0 x 10" 8 . (39) 
Note that although these bounds look more stringent than the corresponding AF = 1 con- 
straints, they scale differently with tan/3 and also involve products of pairs of e^-. Therefore, 
contrary to the AF = 1 case, in principle all of these limits can be evaded for one of the 
couplings by suppressing the other one. Fig. M and Fig. [TU] show the allowed regions for these 
parameters obtained from neutral Higgs contribution to B^ iS —B^ 8 , K-K and D-D mixing 
(see the Feynman diagram in Fig. [8]) . 



9 To obtain a value consistent with the NP analysis of the UTht collaboration, we also used their input for 
computing the matrix elements of the SM AF = 2 Hamiltonian in Eq. ([32]) . 
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B s mixing 



B d mixing 





Re[ef 3 eS*] x 10 1 



-0.6 -0.4 -0.2 0.0 0.2 0.4 



Re[ £ ? 3 €^]x10 10 



FIG. 9: Allowed regions in the complex e^-plane from B^ s -B^ s mixing for tan/3 
m H = 700 GeV (yellow), m H = 500 GeV (red) and m H = 300 GeV (blue). 
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FIG. 10: Allowed regions in the complex ef 2 e2i~pl ane fr° m K—K and D—D mixing for tan/3 = 50 
and m H = 700 GeV (yellow), m H = 500 GeV (red) and m ff = 300 GeV (blue). 



C. Lepton-flavor-violating decays: r — > fj, , r — > e and /i — > e e + e 

In this section, we investigate the constraints that r~ — > [i~fi + n~, t~ — > e~yU + yU~ and 
H —> e~e + e~ place on the flavor changing couplings £39,23^41,13 an d ^21.12; respectively. 
For these decays, the experimental upper limits [83|, |84| are 



B [t~ -)> /i>V] < 2.1 x 10~ 8 
B [t~ -> e~/i + /i^] < 2.7 x 10~ 8 , 
-> e"e + e"] < 1.0 x 10~ 12 



(40) 



at 90% CL. Let us consider the processes r — > n ji + jj, and r — > e fi + fi which are shown 
in Fig. [TTJ The expressions for the branching-ratio for t~ — > e~fi + fi~ can be written as 
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FIG. 11: Feynman diagrams contributing to r~ — > /i~ jjL + /i~ and t~ — > e~ fi + fM~ via neutral Higgs 
exchange. Note that for t~ — > (or fj, — > e~e + e~) two distinct diagrams exist which come 

with a relative minus sign due to the exchange of the two fermion lines. 




300 400 500 600 700 800 300 400 500 600 700 800 300 400 500 600 700 800 



m H [GeV] m H [GeV] m H [GeV] 

FIG. 12: Allowed regions for the absolute value of ef 331 , £2332 an ^ e i2 21 f° r tan/3 = 30 (yellow), 
tan/3 = 40 (red) and tan/3 = 50 (blue) from t~ — > e~[i + fi~, r~ — > and /i~ — > e~e + e~, 

respectively. In each plot only one of the elements or is assumed to be different from zero. 



r — > e 



m" 



tan 4 /3 



12(87r) 3 r r m% 





2 ( 










[\ £ il\^ 


-141 











(41) 



where, T T is the total decay width of the r-lepton. The branching ratios for r~ — > e~e + e~ 
and /i~ — > e~e + e~ can be obtained by an obvious replacement of masses, indices and total 
decays widths. Note that the full expression for general scalar couplings given in Eq. f II 1 6 j) 
of the appendix is different for t~ — > e~fi + fi~ than for r~ — > fi~fi + fi~ and only approaches 
a common expression in the limit of large tan /3 and large Higgs masses. 



Comparing the type-Ill 2HDM expression with experiment we obtain the following con- 
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straints on (assuming e^- = 0) 

|4| 2 +|4| 2 < (2.3X10-) 2 ^ g /500GeV y g[^--ye-e+e- 
1121 ' ' — ' V tan/3/50 / 1.0 x 10~ 12 

|4M4f< (4.2X10-) 2 ( m f/„ 5 °°" eV V B K^fSf 1 ■ < 42 > 

l4| 2 + l4| 2 <(3.7xio- T ^ taii/3/50 ; 21xl0 _ s 

These constraints are also illustrated in Fig. [12] for the experimental limits given in Eq. (140]) . 



tan/3/50 J 2.7x10" 
, 2 / m^/500 GeV\ 4 B [t~ — > fi~fi + fi~ 



VI. LOOP-CONTRIBUTIONS TO FCNC PROCESSES 

We observed in the previous section that all elements efj, ef_- (with i ^ j) and e" 221 
must be extremely small due to the constraints from tree-level neutral Higgs contributions 
to FCNC processes. Furthermore, the constraints on ej, and el get even more stringent if 
both of them are non-zero at the same time due to the bounds from AF = 2 processes. 
Nevertheless, the elements e" 3 23 and e^ 32 are still unconstrained because we have no data 
from neutral current top decays. In addition, also the flavor- conserving elements e{ t are not 
constrained from neutral Higgs contributions to FCNC processes. 

In this section, we study the constraints from Higgs mediated loop contributions to FCNC 
observables. First, in Sec. IVI Al we consider the AF = 2 processes, B s -B s , B d -B d and 
K — K mixing and then examine the constraints on e" 3 2 3 and e 31 32 from b — > s(d)^. Also 
e 22 (e 33 ) can be constrained from these processes due to the relative tan/3 enhancement 
compared to m c (m t ) in the quark-quark-Higgs vertices. In this analysis, we neglect the 
effects of the elements efp which means that we assume the absence of large accidental 
cancellations between different contributions. 

Also AF = processes (electric dipole moments) place relevant constraints on the type- 
Ill 2HDM parameter space, as we will see in Sec IVIFI 

A. B s — B s , Bd — Bd and K — K mixing 

For the charged Higgs contributions to AF = 2 processes we calculated the complete 
set of Wilson Coefficients in a general F^-gauge. The result is given, together with our 
conventions for the Hamiltonian, in the appendix. For the QCD evolution we used the NLO 
running of the Wilson coefficients of Ref . 0, 0] • 

For computing the allowed regions in parameter space we used the same procedure as 
explained in the last section. The results are shown in Fig. [13l [14] and [15] and can be 
summarized as follows: B s -B s (Bd-Bd) mixing gives constraints on e 23 (e" 3 ) which are of 



25 



B s mixing B s mixing B s mixing 




-0.3-0.2-0.1 0.0 0.1 0.2 0.3 -1.0 -0.5 0.0 0.5 1.0 -1.5-1.0-0.5 0.0 0.5 1.0 1.5 

Re[e" 23 ] Re^g] Re[e£ 2 ] 



FIG. 13: Allowed regions in the complex e^-plane from B s mixing for tan /5 = 50 and ran = 
700 GeV (yellow), m H = 500 GeV (red) and m H = 300 GeV (blue). 
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FIG. 14: Allowed regions in the complex e^-plane from B^ mixing for tan (3 = 50 and ran = 
700 GeV (yellow), m H = 500 GeV (red) and m H = 300 GeV (blue). 

the order of 10 _1 (10 -2 ) for our typical values of tan/3 and m#. In addition, B^-B^ mixing 
also constrains to a similar extent as B s -B s mixing. The constraints on €33, €32 and e 31 
are all very weak (of order one). Also Kaon mixing gives comparable bounds on Abs [e^] 
and the bounds on Abs [ej^] are of the order 10 _1 . 

B. Radiative B meson decays: b — > S7 and b — > 

The radiative B decay b — > 57 (b — > dj) imposes stringent constraints on the element 
e 23 ( e i3) while also in this case the constraints on €32 (e^) are very weak due to the light 
charm (up) quark involved (see left diagram in Fig. [T6|) . For these processes both a neutral 
and a charged Higgs contribution occur. Since the flavor off-diagonal elements ef 3 23 and 
e 3i 32 are already stringently constrained from tree-level decays we neglect the neutral Higgs 
contribution here. We give the explicit results for the Higgs contributions to the Wilson 
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FIG. 15: Allowed regions in the complex e^-plane from K — K mixing for tan/3 = 50 and run = 
700 GeV (yellow), win = 500 GeV (red) and mu = 300 GeV (blue). The constraints are practically 
independent of tan/3. 




FIG. 16: Left: Feynman diagram contributing to b — > 57 via a charm-loop containing 632- The 
contribution is suppressed, since the small charm mass enters either form the propagator or from 
the charged Higgs coupling to the charm and strange quark. 

Right: Feynman diagram showing a neutral Higgs box contribution to D — D mixing arising if 
and €32 are simultaneously different from zero. 



coefficients governing b — > s(d) r y in the appendix. 



For B — > Xg'-y, we obtain the constraints on the 2HDM of type III parameters e^- by 
using B[B ->■ X s7 ] from Ref. Q (BABAR) and Ref. 0, (BELLE). Combined and 
extrapolated to a photon energy cut of 1.6 GeV, the HFAG value is 



B [B -+ X.71I2U.6Q.V = ( 3 - 43 ± °' 21 ± °- 07 ) x 10 " 4 



(43) 



In order to estimate the possible size of NP we use the NNLO SM calculation of Ref. 48 1 
(again for a photon energy cut of 1.6 GeV) 



B [B -> X sl ] SM = (3.15 ± 0.23) x 10' 



(44) 
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FIG. 17: Allowed regions for efj from B — > X s ^j, obtained by adding the 2 a experimental error 



and theoretical uncertainty linear for tan/3 = 50 and ran 
(red) and m H = 300 GeV (blue). 
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FIG. 18: Allowed regions for efj from B — > X s ^j, obtained by adding the 2 a experimental error 
and theoretical uncertainty linear for tan/3 = 30 and ran = 700 GeV (yellow), ran = 500 GeV 
(red) and m ff = 300 GeV (blue). 



and calculate the ratio 



R 



b ^ B[B->X sl \r 
oxp ~ B[B^X sl ]\ SM 



(45) 



This leads to a certain range for i?^ S7 . Now, we require that in our leading-order calculation 
the ratio 



b _ sl _B[B^X sl ]\ 



2HDM 



theory 



B[B^X sl f M (46) 
lies within this range. In this way, we obtain the constraints on our model parameters e^- as 
illustrated in Fig. [T7] and Fig. [18] 

The analysis for b — > d'-f is performed in an analogous way. In addition we use here the 
fact that most of the hadronic uncertainities cancel in the CP-averaged branching ratio for 
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B — > Xjy 89|, |90[ . The current experimental value of the BABAR collaboration 9l|, [92| for 
the CP averaged branching ratio reads 



B [B -+ X d7 ]|- p >L6GcV 



1.41 ±0.57) x 10" 



(47) 



Here we take into account a conservative estimate of the uncertainty coming from the ex- 
trapolation in the photon energy cut 93[. For the theory prediction we use the NLL SM 



predictions of the CP-averaged branching ratio B{B — \ ^d7)|s 7 >i.6GeV of Ref. |94j,|95|, which 
was recently updated in Ref. 



931 ] and reads 



B [B X d 



iSM 

l_E 7 >1.6GeV 



;i.54^) x io- 



(48) 



After defining the ratios R^f 1 an d ^theory we continue as i n the case of B [B — > X s ^] in 
order to constrain e" 3 . 

As can be seen from FigJTT] and Fig. [TBI the constraints that B — > X^)^ enforces on 
e 23(i3) are stronger than the ones from B s ^) mixing. Even e 33 can be restricted to a rather 
small range. 

While in the 2HDM of type II b — > 57 enforces a lower limit on the charged Higgs mass 
of 360 GeV [49J this constraint can get weakened in the 2HDM of type III: The off-diagonal 
element e% 3 can lead to a destructive interference with the SM (depending on its phase) 
and thus reduce the 2HDM contribution. Lighter charged Higgs masses are also constrained 
from b — > dq but also this constraint can be avoided by e" 3 . 



C. Neutral Higgs box contributions to D — D mixing 

Nearly all the loop-induced neutral Higgs contributions to FCNC processes can be ne- 
glected because the elements involved are already stringently constrained from tree-level 
processes. However, there is one exception: since the constraints on 32 are particularly 
weak (because of the light charm or up quark entering the loop) this can give a sizable effect 
in D—D mixing via a neutral Higgs box 10 (see Fig. [T6l). As we will use e 31 and e 32 in Sec. IVHI 
for explaining the mentioned deviations from the SM prediction in B — >• tv , B — >■ Dtv and 
B — >■ D*tu it is interesting to ask if all processes can be explained simultaneously with- 
out violating D — D mixing. In principle also charged Higgs contributions to D — D mixing 
arise but we find that they are very small compared to the contributions. The explicit 
expression for the Wilson coefficients can be found in the appendix. 

Fig. [T9l shows the allowed regions in the complex e^e^-plane. The constraints are again 



obtained by using the recent UTFit 82[ analysis for the D D system. 



In principle, one can also get contribution to D — > through H® box and penguin contributions if 

the elements c% 2 and c'li are simultaneously non-zero. However, we observe that they are negligible. 
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FIG. 19: Allowed region in the complex e^e^J-plane obtained from neutral Higgs box contri- 
butions to D — D mixing for tan/3 = 50 and run = 700 GeV (yellow), mu = 500 GeV (red) and 
m H = 300 GeV (blue). 



D. Radiative lepton decays : \i — > cy, r — » e^f and r 



^7 



The bounds on e[ 3 31 and e^ 3 32 from the radiative lepton decays r — > e^f and r — > fi'y 
(using the experimental values given in Table [TT]) turn out to be significantly weaker than 
the ones from r~ — >■ [L r n~ and r~ — > e~/i + /i~. Concerning /i — >• e7 we expect constraints 
which are at least comparable to the ones from pT — > e~e + e~ since /i — >■ e7 does not 
involve the small electron Yukawa coupling entering fi~ — > e~e + e~. In fact, using the new 
MEG results 96[ the constraints from fi — > e7 turn out to be stronger than the ones from 
\i~ — > e~e + e~ (see Fig. [20]) . Note that the constraints from [i~ — > e~e + e~ can be avoided if 
fuCii ~ Tie while the leading contribution to fi — > e'-f vanishes for v u e\ 



~-22 ~ 



Process 


Experimental bounds 


S [r -> /i7] 


< 4.5 x 10~ 8 [97, 9^ 


B [r — )■ e7] 


< 1.1 x 10~ 7 [97] 


B [/i ^ e7] 


< 5.7 x 10~ 13 [96] 



TABLE II: Experimental upper limits on the branching ratios of lepton-flavor violating decays. 

In principle, for /i — > e'-f a simplified expression for the branching ratio in the large tan (3 
limit and v <C m# could also be given. However, due to the large logarithm with a relative 
big prefactor (last term of Eq. (196]) ) this is only a good approximation for very heavy Higgses 
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FIG. 20: Allowed region for e{ 2 (left plot) and e 21 (right plot) from n — > ej for tan/3 = 30 (yellow), 
tan /3 = 40 (red) and tan (3 = 50 (blue) . 

and we therefore use the full expression in our numerical analysis. 

We will return to the radiative lepton decays in Sec. I Villi and correlate them to the decays 
—7- \i~ n + /i~ , t~ — > e~ and /i~ — > e~e + e~ . 



E. B s -)■ 

Setting e?- = only the loop induced charged Higgs contribution to B s — > fi + fi~ (and 
Bd —> exist. This contribution (see Eq. ( |24|) ) gets altered in the presence of non-zero 



elements q-, e.g. e 22 . In the large tan/3 limit, the loop induced result in Eq. fl24j) is modified 
to 



ribs 



C b P s 



m b V t iV ts m. 



m c 22 



tan 2 (3 



log (m 2 H /m^) 



(49) 



2 2M^ r mf^/rn 2 - 1 

The resulting constraints on e 22 from B s — )■ are shown in Fig. [21] and the ones from 

Bd — > are found to be weaker. 



F. Electric dipole moments and anomalous magnetic moments 

1. Charged leptons 

The same diagrams which contribute to the radiative lepton decays for ii ^ if also affect 
the electric dipole moments and the anomalous magnetic moments of leptons for ii = if. 
For this reason we use the same conventions as in Eq. (1931) and express the EDMs of leptons 
in terms of the coefficients c[\ of the magnetic dipole operators 0£r in the following way 
(using that for flavor conserving transitions * = c^*) 

4= 2m ti 1m[<% i ] . (50) 
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FIG. 21: Left: Allowed regions in the complex e^-plane from B s — > for tan/3 = 50 and 

ran = 700 GeV (yellow), ran = 500 GeV (red) and mn = 300 GeV (blue). Right: Allowed regions 
in the e^-mH plane from B s — > for tan/3 = 30 (yellow), tan/3 = 40 (red), tan/3 = 50 (blue) 
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FIG. 22: Left: Feynman diagram contributing to EDMs (for i = f) or LFV decays (for i ^ f) 
involving a neutral-Higgs boson. Right: Feynman diagram contributing to EDMs (for i = f) or 
LFV decays (for i ^ f) involving a charged-Higgs boson. 



In SM there is no contribution to the EDMs of leptons at the one-loop level. This is 
also true in the 2HDM of type II, because the Wilson coefficients are purely real since the 
phases of the PMNS matrix drop out in the charged Higgs contributions after summing 
over the massless neutrinos. However, in a 2HDM of type III, one can have neutral Higgs 
mediated contributions to EDMs. Note that there is no charged Higgs contribution to the 
charged lepton EMDs also in the 2HDM of type III because the Wilson coefficients are 
purely real in this case. Comparing the expression for the EDMs in the 2HDM of type III 
with the experimental upper bounds on d e , and d T (see Table HTT1) . one can constrain the 
parameters ef • (or combination of them) if they are complex. 

We observe that while d e enforces strong constraints on the products Im [e^e^-J and 
Im [ef 2 e2i] (see Fig. [23]), d^ and d T are not capable of placing good constraints on our model 
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EDMs 


141 




d T 


Kl 


Bounds (ecm) 


10.5 x 10~ 28 [22] 


1.9 x 10~ 19 [100] 


G [-2.5, 0.8] x 10~ 17 [101] 


2.9 x 10~ 26 [102] 



TABLE III: Experimental (upper) bounds on electric dipole moments. 
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FIG. 23: Allowed regions in the Im [e^dj-J -mj and Im [ei2 e 2i]~ m # planes from neutral Higgs 
contribution to d e for tan j3 = 50 (blue) , tan (3 = 40 (red) and tan /3 = 30 (yellow) . The constraints 
on Im [e e u ] are not sizable. 



parameters. 

Similarly, following the conventions in Eq. ( 193|) . the anomalous magnetic moments 
(AMMs) can be written in terms of as (e > 0) 

4m? 



Re [<%*] . 



(51) 



The discrepancy b etween experiment and the SM prediction for the muon magnetic moment 



a M = (g - 2)/2 is 



103 



104 



105 



10d Il07[ 



Aa, 



«r - < M 



(3±1) x 10" 



(52) 



In the 2HDM of type II, the sum of the neutral and charged Higgs mediated diagrams gives 
the following contribution to a M (for tan/3 = 50 and m H = 500 GeV): 



a 2HDMII ^ 2 j x 10 _13 



(53) 



which is interfering constructively with the SM. Anyway, it can be seen that the effect is 
orders of magnitude smaller than the actual sensitivity and it even gets smaller for higher 
Higgs masses. 
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FIG. 24: Left: Allowed region in the £23" ^32 pl& n 6 from Aa,^ for real values of 623, 632 and tan/3 — 
50, m H = 700 GeV (yellow), m H = 500 GeV (red) and m H = 300 GeV (blue). Right: Allowed 
region in e^-m/f plane from Aa^ for real values of €22 and tan (3 = 50 (blue) , tan (5 = 40 (red) and 
tan (5 = 30 (yellow) . 



Concerning the 2HDM of type III the discrepancy between experiment and the SM pre- 
diction given in Eq. (152]) could be explained but only with severe fine-tuning. One would 
need to allow for very large values of e l 22 which would not only violate 't Hooft's naturalness 
criterion but also enhance B s — > n + fi~ by orders of magnitude above the experimental limit. 
If one would try to explain the anomaly using 633 and £32 (^12 aR d ^21) one would violate the 
bounds from r~ — >■ (/i~ — > e~e + e~ or /i — > ej) as illustrated in Fig. [241 

In conclusion, neither a type-II nor a type-Ill 2HDM can give a sizable effect in a M and 
both models are not capable of explaining the deviation from the SM. 



2. Electric dipole moment of the neutron 



The neutron electric dipole moment d n can also provide constraints on the parameters 
e\y In the SM, there is no contribution to d n at the 1-loop level since the coefficients are 
real. This is also true in the type-II 2HD M. 



Usin g the theory estimat e of Ref . 



of Refs. 



109 



110 



111 



108 



which is based on the QCD sum-rules calculations 



11211 , the neutron EDM can be written as 



d n = (1 ± 0.5) [lA{d d - 0.25rf u ) + l.le(d 9 d + 0.5^)] 



(54) 



where, d u (do) is the EDM of the up (down) quark and d 9 u ^ define the corresponding 
chromoelectric dipole moments which stem from the chromomagnetic dipole operator 



U R{L) 



m qi q f a^T a P Rmi G% . 



(55) 
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FIG. 25: Allowed regions in the Im [^n 22] _m #" planes from the electric dipole moment of the 
neutron for tan (3 = 50 (blue) , tan (3 = 40 (red) and tan /3 = 30 (yellow) . We observe that d n can 
not provide good constraints on the real parts of 22 . 
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FIG. 26: Allowed regions in the complex 31 -planes from d n for tan/3 = 50 and ran = 700 GeV 
(yellow), mjj = 500 GeV (red) and ran = 300 GeV (blue). We see that the absolute value of 
can only be large if it is aligned to V u b, i.e. ArgfV^] = Argfe^] ± tt which is very important when 
we consider later B — )• rv. 



Similar to EDMs, the (chromo) electric dipole moments of quarks are given as 



2m„. Im 

Hi 



(56) 



Using the upper limit on d n (see Table HTT|) we can constrain some of e^- (for efj = 0) as shown 
in Fig. [25] and [26j These constraints are obtained for the conservative case of assuming a 
prefactor of 0.5 in Eq. (|5~4"|) . The explicit expressions for c q ^ q ^ stemming from neutral and 
charged Higgs contributions to dqf are relegated to the appendix. Note that for the neutron 
EDM we did not include QCD corrections. 



35 



VII. TREE-LEVEL CHARGED CURRENT PROCESSES 



In this section we study the constraints from processes which are mediated in the SM 
by a tree-level W exchange and which receive additional contributions from charged Higgs 
exchange in 2HDMs. We study purely leptonic meson decays, semileptonic meson decays 
and tau lepton decays. Concerning B meson decays we consider B — > tu, B — > Drv 
and B — > D*tv which are, as outlined in the introduction, very interesting in the light 
of the observed deviation from the SM. We consider in addition D( a \ — > tu, Dr s -) — > fiu, 
K(ir) — > eis, K(tt) — > fiu and r — > K(ir)u and look for violation of lepton flavor universality 
via K(ir) — > ev/K{n) — > [iv and r — > K(n)v / K(n) — > [iv. Even though no deviations from 
the SM have been observed in these channels, they put relevant constraints on the parameter 
space of the type-Ill 2HDM. 

For purely leptonic decays of a psudoscalar meson M (and also tau decays to mesons) to 
a lepton £j and a neutrino v (which is not detected) the SM prediction is given by 

2 



B SM [M -+ V] 



m M 



G 2 F m] t m /mIK/^I 



)»7> 



711 



y 1 ^ em 



(57) 



where stands for channel dependent electromagnetic corrections (see Table HV|) . txim 
is the mass of the meson involved and m Uf (m^J refers to the mass of its constituent up 
(down) type quark. The expression for r — > Mv differers by the exchange of the meson 
masses (life time) with the tau masses (life time) and by a factor of 1/2 stemming from spin 
averaging. 

NP via scalar operators can be included very easily: 



B 



NP 



B. 



SM 



III 



M 



c 



Ufdi ,£j 
R 



(m Uf + m di ) m l3 



with 



C 



SM 



4GVK 



2. 



(58) 



(59) 



All quantities in Eq. (1581 are understood to be at the meson scale tum- Like for B s — >■ 
yU. + /i~, the SM Wilson coefficient is renormalization scale independent and the scalar Wilson 
coefficients evolve in the same way as the quark masses. 

In the 2HDM III the Wilson coefficients C% ,at , j and Cg* ' J are given by (neglecting 
terms which are not tan (5 enhanced) 




e kj 



j=i / V k=i 



(60) 



j=l \ k=l 
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Ratio 


Experimental value 


SM prediction 


EM 


B [K -> ei/] /B [K -> H 


(2.488 ± 0.013) x 10" 5 


(2.472 ± 0.001) x 10~ 5 


-0.0378 ± 0.0004 [61] 


B [A" -> H /B [vr -»• H 


(63.55 ±0.11) x 10~ 2 


(63.48 ± 1.37) x 10~ 2 


-0.0070 ± 0.0018 [113] 


B [if -> ev) IB [it -> ei/] 


(1.285 ±0.008) x 10" 1 


(1.270 ±0.027) x 10" 1 


-0.0070 ± 0.0018 [113] 


B [7T -> ezv] /B [7T ->■ 


(1.230 ±0.004) x 10~ 4 


1.234 x 10~ 4 


-3.85% [114] 


B [r -> ifi/] /B [r -»• vrz^] 


(6.46 ±0.10) x 10~ 2 


(6.56 ±0.16) x 10~ 2 


0.0003 ± 0.0044 [1151 


B [r — >■ 7Tf] /B [7T — > /w] 


(10.83 ±0.06) x 10~ 2 


10.87 x 10~ 2 


±1.2% [114] 


B [r -> AV] /B [if -> H 


(1.102 ±0.016) x 10~ 2 


1.11 x 10~ 2 


±2.0% [114] 



TABLE IV: Experimental values, SM predictions and electromagnetic corrections (in the SM) for 
the ratios of charged current processes. The experimental values are obtained by adding the errors 



of the individual branching ratios given in Ref. 7l[ in quadrature. The SM predictions include 
the uncertainties from an d (if involved) as well as the uncertainties due to CKM factors and 
decay constants. As always, we add the theory error linear to the experimental ones. 



Note that C L " 3 is only proportional to one power of tan f3 while C R 1 ' 3 is proportional 
to tan 2 /3. The Hamiltonian governing M —> IjV (r — > Mv) and the Wilson coefficients for 
general scalar interactions are given in the appendix. It is important to keep in mind that, 
since we are dealing with lepton flavour-violating terms, we must sum over the neutrinos in 
the final state because the neutrino is not detected. Note that we did not include the PMNS 
matrix in both Cg 4 ^ 1 ' 3 and ' 3 for simplifying the expressions, since it cancels in the 
final expression after summing over the neutrinos. 

For semileptonic meson decays B — > Dtv and B — > D*tv , which have a three-body 
final state, both the SM prediction and the inclusion of NP is more complicated, as will be 
discussed in subsection IVII A ll 



A. Tauonic charged B meson decays: B — > tv , B — > Dtv and B — > D*tv 

As discussed in the introduction the BABAR collaboration performed an analysis of the 
semileptonic B decays B — > Dtv and B — > D*tv using the full available data set 
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They find for the ratios 



(with 



Hp { * ] ) = B(B D^tu)/B(B ->• D^iu) 

e, \i) the following results: 

Hp) = 0.440 ±0.058 ±0.042, 
Tip*) = 0.332 ±0.024 ±0.018. 



(61) 



(62) 
(63) 



Here the first error is statistical and the second one is systematic. Comparing these mea- 
surements to the SM predictions 

n SM p) = 0.297 ±0.017, (64) 
KsuP*) = 0.252 ±0.003, (65) 

we see that there is a discrepancy of 2.0<r for Hp) and 2.7 a for Hp*) . For the t heory 



predictions we used the update d results o 
based on the results of Refs. 



which rely on the calculations of Refs. 



55 



116) 



1171 . Ilia . Ill9l . Il20l . Il2l| . The measurements of both ratios 



Hp) and H(D*) exceed the SM prediction, and combining them gives a 3.4 a deviation 
from the SM [8|, |9j expectation. 

This evidence for the violation of lepton flavour universality in B — > Dtv and B — > 
D*ru is further supported by the measurement of B — > rv by BABAR 10, |llj and 



BELLE |l2j . Until recently, all measurements of B — > rv (the hadronic tag and the leptonic 
tag both from BABAR and BELLE) were significantly above the SM prediction. However, 
the latest BELLE result for the hadronic tag [3] of B[B tv\ = (0.72±°;^ ± 0.11) x 10^ 4 
is in agreement with the SM prediction [l3 |: 



£>sm [B 



(0.796 



+0.088\ 
-0.087J 



x 10" 



Averaging all measurements, one obtains the branching ratio 

B exp [B -> tv\ = (1.15 ± 0.23) x 10' 



(66) 



(67) 



which now disagrees with the SM prediction by 1.6 a using V u b from the global fit 14]. 

Combining Hp) , H{D*) and B — >• tv , we have evidence for violation of lepton flavor 
universality. Assuming that these deviations from the SM are not statistical fluctuations or 
underestimated theoretical or systematic uncertainti es, i t is interesting to ask which model 
of new p hysics can explain the measured values 
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FIG. 27: Feynman diagram showing a charged Higgs contributing to B — > tv and B — > D^*'tv 
involving the flavour changing parameters and e 32 which affect B — >■ tv and B — > D^*'tv, 
respectively. 

i. 5 Dtz/ and 5 -> L>*tz^ 

Let us first consider the semileptonic decays B — > Dtv and I? — > D*tv . Here the 
Wilson coeffi cients Cr ,t and Cjf' r affect B — >• Drz/ and Z? — >• D*tv in the following way 



54 



55 



132 



11(D) = n SM {D) I 1 + 1.53ft 

n(D*) = k sm (d*) 



^ R "T O r 



+ 1.0 



SM 



c 



cb ,r 



c 



cb ,T 



n 



c6 ,r 
SM 



+ 0.05 


£fCft,T £jCb ,T 


/~icb,T 
°5M 





(68) 
(69) 



For our analysis we add the experimental errors in quadrature and the theoretical uncertainty 
linear on top of this. There are also efficiency corrections to TZ(D) due to the BABAR 
detector jsj which are important in the case of large contributions from the scalar Wilson 
coefficients C c ^[ (i.e. if one wants to expl ain 1Z(D) with destructive interference with the 



SM contribution). As shown in Ref. 1221 ] . these corrections can be effectively taken into 

by an approximate factor of 



account by multiplying the quadratic term in C c ^'[ of Eq. 
1.5 (not included in Eq. ( 16"%]) ). 

Since ef 3 contributes to C^ ,T (the same Wilson coefficient generated in the type-II 2HDM) 
it cannot simultaneously explain TZ(D) and 1Z(D*) . Therefore, we are left with e 32 , which 
contributes to B — > Dtv and B — > D*tv . In the left frame of Fig. EH] we see the allowed 
region in the complex e^-plane, which gives the correct values for TZ(D) and 1Z(D*) within 
the 1 u uncertainties for tan f3 = 50 and nin = 500 GeV, and the middle and the right 
frames correspond to the allowed regions on e 31 from B — > tv . 
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FIG. 28: Left: Allowed regions in the complex e^— plane from TZ(D) (blue) and 1Z(D*) (yellow) for 
tan /3 = 50 and mu = 500 GeV. Middle: Allowed regions in the complex e^-plane combining the 
constraints from B — )■ rv (1 a (yellow) and 2 a (blue)) and neutron EDM (green) for tan/3 = 50 
and vrtu = 500 GeV. Right: Allowed regions in the m#— plane from B — >■ rv for real values of 
e^i and tan (3 = 50 (green) , tan (3 = 30 (orange) . 
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FIG. 29: Left: Allowed regions in the mjj- plane from B — >■ rv for real values of q 3 and 
tan/3 = 30 (yellow), tan/3 = 40 (red), tan/3 = 50 (blue). Right: Allowed regions in the complex 



e 13' e 23 



and e^-planes from B 



rv for mH = 700 GeV (yellow), mjj = 500 GeV (red) and 



mi? = 300 GeV (blue). 



2. B^rv 



In principle, B — >■ rv can be explained either by using (as in 2HDMs with MFV) or by 
(or by a combination of both of them). However, €33 alone cannot explain the deviation 
from the SM without fine tuning, while is capable of doing this 16]. 



B —7- rv can also be used to constrain e{ 3 , e| 3 and e| 3 as illustrated in Fig. |29j In order to 
obtain these constraints, we assumed that all other relevant elements (eg 3 and e^, ) are zero. 
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FIG. 30: Left: Allowed region in the m^-e^-plane (for real values of e^) obtained by combining 
the constraints from D — > fj,v and D s — > [iv for tan /3 = 30 (yellow), tan (3 = 40 (red) and tan (3 = 50 
(blue). While the upper bound on e?^ comes from D s — > fj,v, D — > fj,v is more constraining for 
negative values of e^n- The bound on the imaginary part of are very weak. The constraints 
from D s — > tv turn out to be comparable (but a bit weaker) while the ones from D —> tv are 
weak. 



Process 


Experimental value (bound) 


SM prediction 


B [D a -> tv] 


(5.43 ±0.31) x 10~ 2 


(5.36±°J 4 ) x 10~ 2 


B [D s -> H 


(5.90 ±0.33) x 10~ 3 


(5.5Q±g:g) x io- 3 


B [D -> ri/] 


< 1.2 x 10~ 3 


(1.10 ±0.06) x 10~ 3 




(3.82 ±0.33) x 10" 4 


(4.1518J?) x ID" 4 



TABLE V: Experimental values (upper bounds) and SM predictions for Dr s \ — > tv and Di s \ — > fiv 
process es. The SM pr ediction for D s — > \xv mode takes into account the EM correction effects of 



±1.0' 



13c 



134 



139|]. 



B. D 



(s) 



tv and D 



Previously, there were some indications for NP in £) s — > tv 



133 



13J, |l35|. However, 



using the new experimental values for B \D , — > t v] (s ee Table [V} and the improved lattice 
determination for the decay constant fo s 136L Il37| we find agreement between the SM 



predictions and experiment. Nevertheless, it is interesting to consider the constraints on 
the 2HDM of type III parameter space. Charged Higgs contributions to D( s ) — >■ tv and 
-D( s ) — > fiv have been investigated in Ref. 
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FIG. 31: Allowed regions in the m/f-e^ i2 -plane from K — > fiv/ir — > \w and K — > ev/n — > ev 
for real values of % i2 and tan/3 = 30 (yellow), tan/3 = 40 (red) and tan/3 = 50 (blue). The 
constraints are weaker than the ones from K (it) — > ev/K(ir) — > [iv and r — > K(ir)v / 'K(tt) — > \w 
but cannot be avoided assuming the MFV limit (^- = t 1 ). 

The most important constraints on the 2HDM of type III parameter space are the ones 
on (shown in Fig. 1301). Dr s -\ — > rv and Dr^ — > fiv constrains Re [e^] while the constraints 
on Im [e^] are very weak. In principle, also the ratio — > tv/D^ — > fiv could be 
used for constraining deviations from lepton flavor universality, but the constraints from 
K(n) — > ev/ K{n) — > fiv and r — )■ K(ii)v / K(ii) — > fiv turn out to be stronger. 



C. K — > fiu/ir — > iiv and K — > ev/n — > ev 

The ratio RK e2 ,n n = B[K —> £v] /B [tt — > lv\ (£ = e,/i) is useful for constraining e 22 , efi 
and e| 2 because the ratio of the decay constants /#/ 'f n is known more precisely than the 
single decay constants 



61]. 



For obtaining the experimental values we add the errors of the individual branching 
ratios in quadrature and the SM values take into account the electromagnetic correction. 
The corresponding values are given in Table. IIVI The errors are due to the combined un- 
certainties in fx/f-K, the CKM elements and the EM corrections. We obtained the value 
for V us from K — > iriv (which is much less sensitive to charged Higgs contributions than 
K — > fiu/ir — > fiv) and V u d by exploiting CKM unitarity. 

Fig. 1321 illustrates the allowed regions for e 22 by combining the constraints from K — > 
{iv I ix — > fiv and K — > ev/ir — > ev. Like in Di s ^ — > rv and D( s ) — > fiv the constraints are on 
the real part of e 22 while the constraints on the imaginary part are very weak. Concerning 
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FIG. 32: Left: Allowed regions in the m#-e| 3 -plane from r — > Kv/t — > irv. Right: Allowed 
regions in the m^-e^-plane obtained by combining the constraints from K — > [iv/n — > [iv and 
K — > ev/ir — > ev for real values of e 22 . in both plots tan/3 = 30 (yellow), tan/3 = 40 (red) and 
tan (3 = 50 (blue). 

e\y and e\ 2 the constraints from K(tt) — > ev/K(7r) — >■ /i^ will turn out to be more stringent 
but the latter ones can be avoided in the limit — = ^f- (see Fig. [31] and Fig. 133]) . 



D. r — > Kv/t —> nv 

The r is the only lepton which is heavy enough to decay into hadrons. The ratio 
B [t — t- Ku]/B [r — > ttu] can be considered for putting constraints on e 21 , ef 2 and e £ i3 . 

The experimental and theoretical values for this ratio are given in Table. IIV1 We observe 
that the constraints from D° — > fi + fi~ and D — D mixing on e 2l and — > fi + fi~ on ef 2 
are too stringent so that no sizable effects stemming from these elements are possible. Also 
concerning ef 3 , as we will see in the following sections, the constraints from r — > ttu/tt — > fiu 
will be stronger but again the latter ones can be avoided in the MFV limit — — = ^f- (see 

3 33 

Fig.[32]). 

E. Tests for lepton flavour universality: K{ir) — > ev/K(Ti) — > [iv and 
r — > K (tt)v I K (it) — > [iv 

Ki 2 {K y iv) decays (£ = e,fi) are helicity suppressed in the SM and suffers from 
large theoretical uncertainties due to the decay constants. However, considering the ratio 
Rk 12 — B[K — > eu] /B[K — > jiu] the dependence on decay constants drops out. 
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FIG. 33: Allowed regions in the mij—efj— plane obtained by combining the constraints from K — > 
ev/K — > fiv and ir — > ev/ir — > fiv for real values of e\- and tan/3 = 30 (yellow), tan/3 = 40 (red) 
and tan/3 = 50 (blue). The constraints on (affecting the electron coupling) are more stringent 
than the constraints on e\ 2 (which affect the muon coupling). 

In the 2HDM of type II the charged Higgs contributions to K(ir) — > ev/K(n) — > fiv and 
r — > K (ir)v/K (7r) — > fiu drop out. This is also true in the 2HDM of type III (for q,- = 
with i 7^ j) as long as the MFV-like relation e^/m^tli/rne is not violated. 



1. K — > ev/K — > fiv and ir — > ev/ir — > /iv 



K — > ev/K — )■ fiv is a very precise test of lepton flavor universality 1401 ] (see table. II VI) . 
Including NP entering via scalar operators modifies this ratio according to Eq. (158|) . 

We find strong constraints on e\ 2 (which affect the coupling to the muon) and the con- 
straints on e\ x (where the coupling of the electron is involved) are even more stringent. Like 
for D( s ) — > rv and D( s ) — > fiv the constraints are much better for the real part of than the 
imaginary part. Note that these constraints are obtained assuming that only one element 
efj is non-zero. In the case e e 22 / ' m n =e iil ' m e where lepton flavor universality is restored no 
constraints can be obtained. 

Alternatively, the ratio ir — > evjix — > fiv can test lepton flavor universality. We find that 
the constraints from ir — > ev/n — > fiv are comparable with the ones from K — > ev/K — > fiv. 
Our results are illustrated in Fig. 
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FIG. 34: Allowed regions in the m#-q 3 -plane from r — > -ku/tt — > \iv for real values of e i3 and 
tan/3 = 30 (yellow), tan/3 = 40 (red), tan/3 = 50 (blue). The bounds on the imaginary parts are 
very weak. 

2. t — > Ku / K —7- \iv and r — > ixv j-K — > [iv 

The ratios r — > KvjK — > fiu and r — > ttv/tt — > fiu are very similar to K(tt) — > 
ev/K{jx) — > fiu: all dependences on decay constants and CKM elements drop out and 
they are only sensitive to NP which violates lepton-flavour universality. The corresponding 
experimental and the theoretical values for these ratios are given in Table. IIVI 

We find that the constraints on e e i3 from r — >■ tiu/tt — > fiu are stronger than the ones from 
t — > Kv / K — >■ fiv and they are shown in Fig. [34j 



VIII. UPPER LIMITS AND CORRELATION FOR LFV PROCESSES 

In Sec. [V] we found that the neutral current lepton decays r~ — > fi~ fi + fi~ and r~ — > 
e~fi + fi~ give more stringent bounds on the elements e e 32 23 an d 13 than the radiative 
decays r — > /I'j and r — > ej. Also the LFV neutral meson decays B s ^ — > rfi, B s ^ — > re, 
B s ,d A*e cannot be arbitrarily large in the type-III 2HDM due to the constraints from 
B s ,d and r~ — > r~ — > e~fi + fi~, fi~ — > e~e + e~ (assuming again the 

absence of large cancellations) 11 . 

Therefore, in this section we study the upper limits on B s ^ — > r/i, B s ^ — >■ re, B s ^ — > fie 
and the correlation among — > r~ — > e~/i + /i~, ji~ —> e~e + e~ and r — > ^7, 

r — > ej, /1 — > e-f in the type-III 2HDM. 



11 



see e.g. Ref. 14ll 1 14=21 Il43j for an analysis of NP in B s ,d —> Tfi. 
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A. Neutral meson decays: B s d — > t/j,, B s 



re and B s d — >■ fie 



In the SM (with massless neutrinos) the branching ratios for these decays vanish. Also in 
the 2HDM of type II these decays are not possible (even beyond tree- level). In the type-Ill 
2HDM, these decay modes are generated in the presence of flavor-violating terms efj and 
there exists even a tree- level neutral Higgs contribution to B s — > £f£j (Bj — > £f£j) if also 



-23,32 



T^O(65WO). 



Observables 


B(B S -> fie) 


B(B d -> fie) 


£(£? d -> r/i) 


-> re) 


Upper bounds 


2.0 x 10~ 7 [144] 


6.4 x 10" 8 [144] 


2.2 x 10~ 5 [145] 


2.8 x 10~ 5 [145] 



TABLE VI: Upper limits (90 % CL) on the branching ratios of the lepton flavor-violating B meson 
decays. 

In the large tan/3 limit, v <C m# and neglecting the smaller lepton mass, the correspond- 
ing expressions for these branching ratios take the simple form 



B [B q iflj] « JVj 
with q — d,s, = Nf- and 



tan 0/50 \ ' 
m^/500 GeV J 



-</3| 



I I 3g | 



(70) 



TV?! « 2.1 x 10 7 — , 

21 0.229 GeV ' 

N& « 1.6 x 10 7 ^ , 

21 0.196 GeV 



A" o 2 « 1.7 x 10 7 ^ — 

31,32 0.229 GeV 



Nt o 2 « 1.2 x 10 7 ^ — 

31,32 0.196 GeV 



(71) 



Note that the expressions for the branching ratios are not symmetric in and e^. Since 
experimentally both B q — > £f£j and B q — >■ are combined we compute the average 



(B [fi, ^7] + B [B 



])/2 



In order to obtain the upper limits we insert the biggest allowed values for Abs 
For 32 (ef 3 31 ) we use the biggest allowed absolute value compatible with the bounds from 
B s -> fi + fi~ (B d -»■ fi + fi~). As we can see from Fig. U] (Fig. [5]) the absolute value for £32 
(e^) can be bigger than (ef 3 ). For the leptonic parameters ef 331 and e2 3 32 we use the 
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S[Bs^T/2] 





FIG. 35: Upper limits on the branching ratios of the lepton flavor violating B meson decays as a 
function of run for tan /3 = 30 (yellow) , tan (3 = 40 (red) and tan /3 = 50 (blue) . 

constraints obtained from r~ — > fi~fi + fi~, r~ — > e~/i + /i~ (see Sec. IV CP 

V I < 42xl0- 3 r m ^ 50 ° GeV V 
|e slils | < 4.2x10 ^ tan/?/5o j , 

|4 23 | < 3.7x10-3 MOGeVy 

1 32 ' 231 - V tan/3/50 J ' V ; 

while for ef 2 21 we use the combined constraints from /i~ — > e~e + e~ and from /1 — >■ e7 (see 
Sec. lVTD]) . 

Our results are shown in Fig. [35j We see that for bigger Higgs masses larger values for 
the branching ratios are possible. 



B. Radiative lepton decays: r — > ^7, r — > ej and [i — > ej. 

In Sec. IV CI and Sec. IVIDI we found that the radiative lepton decays r — > ji'y and r — >■ 
give less stringent bounds on the parameters 63332 and ef 331 than the processes t~ — > 
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FIG. 36: Left: 
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g ^e-e+e- l as a f unc ti° n of m# assuming that only e\ 2 (red) or e 2 \ (blue) is 



different from zero for tan /3 = 50. Right: gpr^=^^p] as a function of m# assuming that only 

(red) or e^ 2 (blue) is different from zero for tan f3 = 50. 
For scenarios in which both and e| 2 (^12 an d e|i) are different from zero the 2HDM of type III 



predicts the ratio 



g[r-»/£r] 



to be within the yellow region. These ratios are 

£>[t— >erf\ 



B[t~— \ B[ti- —>-e-e+e- 

to a good approximation independent of tan /3 for tan f3 > 20. The behavior of 
shown here) is very similar to the case of 3 — > 2 transitions. 



e-^j (not 



\i H + H and r — > e /i + fi while the constraints on e[ 2 21 from ii — > are stronger than 
the ones from ii~ — > e~e + e~ . 

There are however interesting correlations between these decays in the type-Ill 2HDM. 
In the large tan j3 limit and for v <C rrin we obtain the following relation 



'Si i 



] 





- 4 


HI 


<7fl 


' + 4 


\ € fi 


IJ 


24tt 




2 ( 


\ e if 


|2 1 
+ 


e l 1 





(73) 



As already noted in Sec. IVI D\ we stress that this formula is only a good approximation 
for very heavy Higgs due to the large logarithmic term in the expression for li — y iff (see 
Eq. (19(51) ). Therefore, the relation in Eq. (173"]) gets modified for lighter Higgs masses as shown 
in Fig. [33 We see that, as expected, for very large Higgs masses the ratios approach 



B [4-^/7] 



a Trie 

"em t 



7 i i 



24vr m|. 



B [4-^/7] 



"em « 



'S i 3 



67r mf . 



for 4 f ? 



(74) 



for e% ± 



where, we assumed that ej^/efi = m ijl m ii an d that only one flavor changing element e^, % 
is different from zero. 
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IX. CONCLUSIONS 



In this article we studied in detail the flavor phenomenology of a 2HDM with general 
Yukawa couplings. Motivated by the fact that the 2HDM of type III is the decoupling limit 
of the MSSM we assumed a MSSM-like Higgs potential. In our analysis we proceeded in 
several steps: 

1. We gave order of magnitude constraints on the parameters efj from 't Hooft's natu- 
ralness criterion and found that all couplings except e" 3 3i and e\\ 2 2 should be much 
smaller than one. 

2. Considering tree-level FCNC processes we constrained the elements ef, (i ^ j) and 
e i2 21 fro m neutral meson decays to muons and from AF = 2 processes, finding that 
they are tiny for the values of m# and tan/3 under investigation (assuming efj = 0). 
In the lepton sector the absolute values of all flavor off-diagonal elements were 
constrained from r~ —> fi~fi + fi', t~ — > e~fi + fi~ and fi~ —> e~e + e~ to be very small. 

3. After having found that the off-diagonal elements must be very small due to con- 
straints from tree-level contributions to FCNC processes we considered charged Higgs 
contributions to K—K , B s —B s , B^—B^ mixing and b — > s(cf)7 arising at the one-loop 
level. In these contributions the so far unconstrained elements e" 33i (and also e^) 
enter for the first time and we found that, setting efj = (with i ^ j), e" 32 3 should 
be rather small. Furthermore, the electric dipole moment of the neutron and of the 
charged leptons constrain e" 1; e^, ^21 > e ti an d e ij; respectively. Respecting all other 
constraints, no sizable effect in a M is possible. 

4. Keeping in mind the constraints from the previous steps, we considered the possible 
effects in charged current processes. Here we found that tests for lepton flavor univer- 
sality constrain the differences e^/m^ — ej^/m^.. Most importantly, the unconstrained 
elements e 31 and e 32 enter the processes B —¥ tv and B — > D^tv directly (with- 
out CKM suppression) and can remove the tension between experiment and theory 
prediction observed in the SM simultaneously. 

5. Finally we gave upper limits on the lepton flavor violating neutral B meson decays in 
the 2HDM of type III and correlated the radiative lepton decays to t~ — > fi~fi + fi~, 
t~ — > e~ and /i~ — > e~e + e~ . 

In Table |VlT] and IVIlT] we list all processes which have been under consideration and quote 
the constraints placed on the parameters eff for our benchmark point m# = 500 GeV and 
tan = 50. 

In summary, combining the constraints from Table IVHI and IVIIII the following bounds 
on the absolute values of the parameters e?- and ef • (for our benchmark point with = 
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Observable 


Results 


Neutral meson decays to muons 


B s ->■ fi+fj," 




ei, 2 < 3.0 x 10~ 5 , 


e^ 3 < 1.9 x 10" 5 , e e 22 < 2.0 x 10~ 3 


B d ->■ 


< 1.1 x 10~ 5 , ef 3 < 9.4 x 10~ 6 




4i < 1.6 x 10" 6 , ef 2 < 1.6 x 10~ 6 


D° -»• 


|^| < 3.0 x 10~ 2 , \ef 2 \ < 3.0 x 10~ 2 


AF = 2 processes 


B s — B s mixing 




4A% < 9-2 x 10- 10 , |e&| < 0.18, |e&| < 1.7, |e&| < 0.7 


B d -B d mixing 




ef 3 eli < 3.9 x HT 11 , |e&| < 0.2, | e ? 3 | < 0.04, |e&| < 1.9 


K — K mixing 


e d l2 4\ < 1.0 x 10~ 12 , |e&| < 0.25, |e| 3 | < 0.14 


D — D mixing 


|ejfee&| < 2.0 x 10" 8 , |e&e&| < 0.02 


Radiative B decays 


b — > sj 


|e&| < 0.024, |e&| < 0.55 




|e? 3 | < 7.0 x 10~ 3 


Radiative lepton decays 


fi — y 67 


e 12 


< 1.7 x 10~ 4 , 4i 


< 2.2 x 10- 4 , 55 < g f^j < 86 

— ' — o[/i — >e e+e \ — 


r — )• e7 


0.19 < m B[ ^ ei \ =, < 0.35 


r — )• /i7 


0.19 < m B[ :^i < 0.35 


Neural current lepton decays 


/x _ — >■ e~e + e~ 




J 

-12,21 


< 2.3 X 10~ 3 


r~ — >■ e~/i + /i _ 




e 13,31 


< 4.2 x 10~ 3 


T _ — )■ [L~ [l^ [l~ 




J 

e 23,32 


< 3.7 x 10~ 3 



TABLE VII: Results obtained in the type-Ill 2HDM from various processes for tan/3 = 50 and 
m H = 500 GeV. 
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Observable 


Results 


Charged current processes 


JD —f TV 


2.7 x 10~ 3 < |e&| < 2.0 x 10~ 2 , 4 


< 6.0 x 10~ 2 


JD — r UTV cv I j — ? 1J TV 


0.43 < |ef 2 | < 0.74 


J-J g 7 1 V 66 J-S(s) * 


|Re[e&]| <0.2 


D TV 




K — > fi{e)v /ir — > fi[e)v 


Re [e* 2 ] 


< 1.0 x 10~ 3 


K(tt) — > ev/K(ir) — > fiv 




Re [4] 


< 2.0 x 10~ 6 , 


Re [4] < 5 -° x 10 " 4 


t — > K(tt)v/K(tt) — > flV 


-4.0 x 10~ 2 < Re [4] < 2.0 x 10" 2 


t — > Kv jr — > TTV 


S3 


< 0.14 


EDMs and anomalous magnetic moments 


d e 


Im 


e 12 e 21. 


< 2.5 x 10~ 8 , 


Im [ e 13 e 31. 


< 2.5 x 10~ 9 






d T 




d n 


(Im^JI < 2.2 x 10~ 2 , |Im [e^ 2 ]| < 1.1 x 10 _1 , Arg[e^] = Arg[V ub ] ± tt 




Deviation from the SM cannot be explained 


LVF B meson decays 


B s -)• r/x 


B [B 8 Tfi] < 2.0 x 10~ 7 


B s -> /ie 


B [B s fie) < 9.2 x 10" 10 


5 S ->■ re 


B [B s re] < 2.8 x 10~ 7 




B[B d Tfi] < 2.1 x 10~ 8 


-Bd ->■ fJ-e 


B [B d fie] < 9.2 x lO" 11 


B d -» re 


B [B d re] < 2.8 x 10" 8 



TABLE VIII: Results obtained in the type-Ill 2HDM from various processes for tan /3 = 50 and 
m H = 500 GeV. 
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500 GeV and tan (3 = 50) are obtained: 



/ 



< 



< 



3.4 x 10" 4 3.0 x 10~ 2 7.0 x 10~ 3 

3.0 x 10~ 2 1.4 x lO" 1 2.4 x 10~ 2 

2.0 x 10- 2 7.4 x lO" 1 5.5 x 10" 1 

1.3 x 10" 4 1.6 x 10~ 6 9.4 x 10~ 6 



1.6 x 10 
1.1 x 10 



-5 



2.6 x 10~ 4 2.0 x 10~ 5 
3.0 x 10~ 5 1.4 x 10" 2 



(75) 



/ 



ij 



< 



2.9 x 10" 6 1.7 x 10" 4 4.2 x 10~ 3 
2.2 x 10~ 4 6.1 x 10" 4 3.7 x 10~ 3 
4.2 x 10~ 3 3.7 x 10~ 3 1.0 x 10^ 2 



These bounds hold in the absence of large cancellations between different contributions. 
Note that in Eq. (175|) we applied the naturalness bounds in case they were stronger than 
the experimental limits. 

It is interesting that B — > tv , B — > Dtv and B — > D*tv can be explained simultaneously 
in the 2HDM of type III without violating bounds from other observables and without 
significant fine-tuning. It remains to be seen if these tensions with the SM remain when 
updated experimental results and improved theory predictions will be available in the future. 
In order to further test the model and constrain the parameters €32 (e^) we propose to study 
H°, A -> tc (H°, A -> tu) at the LHC. 

Acknowledgments : This work is supported by the Swiss National Science Foundation. 
We thank Jernej Kamenik and Ulrich Nierste for collaboration in the early stages of this 
work. A.C. also thanks Jernej Kamenik for useful comments on the manuscript. 



52 



X. APPENDIX 



In this appendix, we collect the Wilson coefficients (to the relevant precision at the 
matching scale) which are needed for the calculation of b — > s(d) , ~f, AF = 2 processes 
(i.e. neutral meson mixing), leptonic neutral meson decays (AF = 1 processes), B — > tv , 
B Dtv , B D*tv , £>( s ) -> £u£, K(ix) ->■ r -> K(tx)v, LFV radiative lepton 
transitions, EDMs of charged leptons and neutron, as well as the AMM of the muon. In 
addition, we give general expressions for some branching ratios, the explicit form of the loop 
functions entering our results and summarize the input parameters used in our analysis in 
tabular form. 



A. Loop functions 

We give the explicit form of the loop functions entering our results. In the limit of 
vanishing external momentum the one and two-point functions are defined as 



1 



(2n) d (k 



d (b2 _ m 2\ ' 



. 2 2 . 167T 2 r d d k 
S (mf,m|) = — : — /i 4 a J 



(76) 



(2n) d (k 2 - mf) {k 2 - ml) 



2\ ' 



where /x is the renormalization scale. 

The loop functions Cq (three-point) and D (four-point) are defined in analogy to B , 
but with three and four propagators, respectively. Evaluating these loop functions yields 
(with d = 4 - 2e) 



A (m 2 



Bo (ml, ml) 



1 H 7s + In (47r) + In 

1 m?hl (^ 

- - lE + In (4vr) + 

e mi 



0(e 
ml In ( At 



m% 



+ 0(e) 



(77) 



2 2 2 \ B (ml, ml) - Bo {ml, ml) 
C (ml,ml,mi) = 



mS — m% 



(m 2 \ ( Tfir\ ( m 2 \ 

— \ + mkn In I — I ) + m\ml In I — | 
m 2j \ m iJ \™>iJ 



[m. 



m% 



m\ 



m 



Do (ml,ml,ml,ml) 



Cq (ml, ml, ml) — C (ml, m\, mf) 



m% 



m- A 



(78) 



Here, the one and the two-point loop functions A , B are UV-divergent and e is the UV- 
regulator. 
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At various places also the functions C 2 and D 2 appear, which have, compared to Co and 
D , an additional factor k 2 in the numerator of the integrand. These functions read 



C 2 (ml, ml, ml) = B (m\,m 2 2 ) + m 2 z C Q (m\,m 2 2 ,ml) , 

D 2 (m\, ml, ml, mf) = C (m\, ml, mf) + ml D (m\, ml, ml, ml) . (79) 
B. Radiative b — > s(d)j decays 

Concering new physics contributions to b — > s(d)j, we work in leading logarithmic (LL) 
precision in this paper. As mentioned before, we use these processes to constrain certain 
elements e^-. For this purpose, we put the e^— couplings (which are already constrained to 
be very small) to zero. When also neglecting the mass of the strange quark and further 
neglecting operators with mass dimension higher than six, we obtain the same effective 



Hamiltonian as in the SM, reading for b — > 57 (see e.g. Ref. 47j ) 



H eff ~ 



^VnVZY^CiOi. (80) 



For b — > d'j the CKM structure is slightly more complicated (see e.g. Ref. 94J). In our 
approximation only the Wilson coefficients C 7 and C 8 of the operators 

7 = -^m b sa^P R bF flu ; 8 = -^m b sa^T a P R bG« v (81) 

get new physics contributions. They are induced through charged Higgs bosons propagating 
in the loop (neutral Higgs boson exchange leads to power suppressed contributions which 
we neglect). For b — >■ 57 the new physics contributions read (with i/j = m 2 l ./m 2 H+ and 
At = V tb V*) 

Np ^J_V^ T RLH ± i, v LRH ±C T,XYiyj) 
7 ~ \ m 2-^ u J d 2 u J d 3 



A* m b z — f 3 3 m u 

j=i 



V 2 r RLH±±-rRLH± C 7,YY(Uj) (RC) -, 



H ■ 

J = l 



n NP _ V<2 1 -pR.LH ± -k -pLRH ± 

ITJ7 ^— » J 



^yy (i/j) 
A t ^ ' ' u > di ml 



v 



EyRLH ± * yRLH ± 
U] d2 Ujd 3 m 2 
j=l U J 



while for b — > d'j the label d 2 and At = V tb V* have to be replaced by d\ and X t = V tb V^, 



respectively. The functions C® XY , C® YY , C\ XY and C\ YY were introduced in Ref. |47| : 
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their explicit form reads 



O® XY^Vj) 



Vi_ 
12 

Cg^XYiVj) — ¥ 

Ci,YY{yj) — ^ 

Cs,YY(Uj) = 



-5y] + 8yj - 3 + (6y 3 - - 4) lny, 

(l/i - I) 3 
-y? + 4y J --3-21ny j 

(Vi " i) 3 

-8y? + 3y] + I2y 3 - 7 + (18y? - 12^) lny, 

(yj - 1) 4 

-I/ 3 + %f ~ jfyj - 2 - 6y 3 In y 3 - 
" I) 4 



(83) 



In Eq. ( )82|) we retained the contributions from internal up- and charm-quarks, although 
these contributions are subleading. 



C. Wilson coefficients for AF = 2 processes 



The extended Higgs sector of our 2HDM of type-Ill also leads to extra contributions to 
AF = 2 processes (B s , Bj, Kaon and D mixing) which can be matched onto the effective 
Hamiltonian 

k%T 2 = E c i °i + E G 'i °i + h - c - > ( 84 ) 

j=l 3=1 

where the operators read in the case of B s mixing 

Ox = (s a ^P L b a ) (sp-fP L bp) , 2 = (s a P L b a ) {s p P L bp) , 

3 = (s a P L b p ) (spPM , 4 = (s a P L b a ) (s p P R bp) , (85) 

C 5 = {s a P L b p ) {spP R b a ) . 

a and /3 are color indices and the primed operators can be obtained from Oi.2,3 by inter- 
changing L and R. Similarly, the corresponding operator bases for Bd, Kaon and D mixing 
follow from Eq. (185]) through simple adjustment of the indices. 

In the following subsections we present the contributions to these Wilson coefficients 
arising from: 1.) one- loop box diagrams with charged Higgs boson exchange; 2.) tree- level 
contributions induced by neutral Higgs boson exchange; 3.) box diagrams involving neutral 
Higgs bosons, relevant in the case of D mixing. 



1. Charged Higgs box contributions 



For definiteness, let us consider B s mixing. The corresponding Wilson coefficients for 
Bd and Kaon mixing follow by a simple adjustment of the indices. We have performed our 
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calculation in a general gauge. The non-vanishing Wilson coefficients from pure charged 
Higgs boxes are given by 



,m H± , , 



up " b u k -i "°H±i m H± j ' 



1 3 

r< _ vR L H ± *-pR L H ± -pRL H^-pRL H± D ( m 2 m 2 m 2 

~ 128vr 2 • _ Ujd2 Ujds Ukd2 Ukd3 \ up Ufc ' H± 

1 3 

n — ~ L V m m T LRH± *T RL H± T LRH± *T RL H± D „ ( m 2 m 2 m 

~ 32tt2 ?-L 3 k u > d2 Ujds Ukd2 Ukd3 { m up m u^ m 

1 3 

C 4 = — 5 £ m u .m Ufc r^ H± *r^ H± r^/ ± *rf'^ ± J Do fm?, ,m 2 ,m 

4 i c_2 ' U J n fc "j<Z2 "ja3 «fc<Z2 «fea3 u V ' «fe' 

±u/ ' j,k=l 

1 3 

/nr _ -pLR H ± *-pLR H* 1 -pRL H ± *-pRL H ± pi 1^2 _,2 ™2 ™2 "\ 

~~ 327T 2 " • _ " jd2 Ufed2 " fed3 V U i ' ' H± ' H± J ' 

j.k — 1 



H± ,m H± , , 



3 

s 

0,' 



The sum of the charged Higgs— W and charged Higgs— Goldstone-boson boxes is given by 

,2 3 



(87) 



r — ^ 2 ( m <m T/*T/ V RLH± V RLH± * 

4M 2 , D (ikf^, m^ ± , , m* fc ) - D 2 (ikf 2 ,, m 2 ^ , m 2 ^ , r<) 
X 4M| 

1 2 3 

r» _ ^2 ft/ T/* 'pLRH ± *-pLRH ± 

° 4 " le^T^A J fc2 ^ 

X 

We stress here, that we want to use B s mixing only to constrain certain e^— couplings, 
because the e^— quantities are already contrained to be very small. We therefore only 
took systematically into account those contributions to the Wilson coefficients which stay 
different from zero in the limit efj — > 0. At first sight, the Wilson coefficient C4 seems to 
be gauge dependent. However, when using the unitarity of the CKM matrix (entering the 
expression for C4 both, explicitly and implicitly through the T— quantities), we find that the 

lich we put to zero in our 



146] . The only difference 



£— dependent terms are always proportional to an element efp w 
analysis. Also note that our result agrees with the one of Ref. 
is that we neglected gauge dependent terms corresponding to higher dimensional operators. 
The Wilson coefficients of the primed operators can be obtained by interchanging L and R 
in the corresponding unprimed ones. 
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2. Tree-level H2 contribution 



The Wilson coefficients from neutral Higgs mediated tree-level contributions to B s mixing 
read: 



a 



2 ~ 9rr,2 ^ d 2 d 3 



r H k - ST - 1 r LRH° LRH°* 
fe=l ^ 

The corresponding coefficients for Bj, Kaon and D mixing follow by a careful adjustment 
of the indices. 

Note that in the limit of large tan (3 and 3> v , C 2 k and C 2 k vanish and we only get 

H° 

a contribution to C 4 k . 



3. Neutral Higgs box contribution to D mixing 

The Wilson coefficients resulting from the neutral Higgs box contribution to D mixing 
are given as 

3 3 



— 1 ^ ^ LRHk!* LRH l 2 LRH k 2 * LRH k 1 ( 2 2 
= 1287T 2 ^ ^ ^ U2 "« 1 2 r« 2 «i 2 2 r«i«i 2 1 £>2 (j 71 ^ 5 m Uj 

31 J2 = l fel,fc 2 =l 

3 3 

— 1 ^-^ ^-^ LRH^-k LRH° 2 * LRHl * LR H% * 

^ 2 — tttt-^ 2^ ^u^^^ru^ui 1 T U2Ujl 2 T Uj2U1 1 T U2Uj2 2 



2 2 



32vr 2 

iu'2=i fci,fc 2 =i 



C 3 = 0, 

3 3 



— 1 ^-^ ^-^ LRHl^k LRH° 2 * LRH° 2 LR i?° 



167T 

jl,j 2 =l fcl,fc 2 =l 



x (m 2 ,m 2 m 2 H o ,m 2 H? ) , (89) 



3 3 

— 1 ^ ^ LRHk!* LRH l 2 LRH l 1 LRH k 2 * ( 2 2 

= 1287T 2 ^ ^ 1 ^ il ii«2 2 r ui « i2 1 r n2nj2 2 D2ym u .^,m u 



128tT 2 

jl J2 = l fci,fc 2 = l 



2 2 
, 777. j^o , ?77^o 

1 fei k 2 



The indices ji,j2 describe the internal up- type quarks while ki,k2 stand for neutral Higgs 
indices (H°,h°,A°). Moreover, the primed Wilson coefficients can be obtained from above 
by the replacement L -H- R in the couplings. 
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D. Semileptonic and leptonic meson decays and tau decays: B — > (D^)tv, Dr s \ — >■ 
lug,, K(-k) — > tvi and r — > K(tt)v processes 

These processes are governed by the effective Hamiltonian 

g. £j u fdi,ij Q u fdi/j £,Ufdi,£j QUfdi,£j ^ ^jUfdi,lj ^Ufdi,lj _|_ q (90) 

with the operators defined as 

^ = Uf p Rd .l.p LVi (91) 
0"&*i =u f P L d i £ j P L u. 

Here, for tauonic B meson decays £j = r, di = b and Uf = u (v,f = c) for B — > tv [B — > 
Dtu and B — > D*tu ). For D s — > IjV (D — > £ji>), u f = c and di = s (d), for r — > K(tx)u, 
£j = r, Uf = u and di = s (d) and for K(n) — > £jV we have £j — e, Uf — u and di = s (d). 
The Wilson coefficients in 2HDM of type III at the matching scale read 



°SM = 


AG F 
V2 


^Ufdi i 




R. 


-1 


yLR^ 

Ufdi 


-pLRH^-k 


m 2 H± 




-1 


-pRL H^ 1 

Ufdi 


-pLRH^-k 


m 2 H± 



(92) 



E. Lepton flavour violation (LFV): ii — > £fy processes 

The radiative lepton decays £i — > ifj (£ — e, \i or r) are induced by one- loop penguin 
diagrams with internal neutral or charged Higgs bosons. The result for the one-loop decay 
amplitude can be written as a tree-level matrix element of the effective Hamiltonian 

Keff = c l t d£ + cT O't , (93) 

where and c/ £ ' are the effective Wilson coefficients of the magnetic dipole operators 

= rnzJ f cT»„P R{ L)£iF^ . (94) 

With these conventions, the branching ratio for the radiative lepton decays £{ — > iff reads 



,5 



BV^ia\ = £^(\4 li \ 2 + \cl"\ 2 )- 
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The neutral Higgs (H® = H°, h°, A ) penguin contribution to c^" 1 is given by 



C RH° 



3 

E 



^LRH° LRH°* m £ f v LRHl* v LRHl 
" ' J- f.f . r 1 « .« *- p .f. 
3 TTl£ 



m tj LRH°LRH° l n , / m ^ 

- r Vj r ^ I 9 + 61n 



™4 



??2 



(96) 



and c^ 1 can be obtained from c^ 1 by interchanging L and i?. Similarly, for the charged 
Higgs penguin contributions we find 



3847r 2 m 2 



H± 3=1 



2^ 1 Vj ti 1 Vj i f 



f c T^LRH ± -pLRH ± * 

R H± ~ m ti 3847r 2 m^ ± ^ ^ v & 



(97) 



F. Wilson coefficients for EDMs and the anomalous magnetic moment of the muon 

1. Wilson coefficients for EDMs of charged leptons and the anomalous magnetic moment of 
the muon 



As in the case of the LFV processes discussed in the previous section, we again have 
both neutral and charged Higgs penguin contributions to the flavor conserving radiative 
transitions l{ — > Icj. The corresponding effective Hamiltonian is obtained from Eq. ( )93l) 
and Eq. (194]) by identifying 1$ with 1^. The contribution to the effective Wilson coefficients 
related to neutral Higgs bosons (propagating in the loop) reads 



R H ° 



k,j=i 



1927r 2 m 2 



m tj LRH°LRH° I [ m i 3 

-IY„. T„.„. I 9 + 61n 1 



rn 



H° h . 



Mi* 



'L H° — U R H°i 



while the charged Higgs penguin contribution leads to the (real) coefficients 

3 



C L H± ~ C R H± 



yLRH ± yLRH ± * 



(98) 



(99) 



(100) 



2. Wilson coefficients for neutron EDM 



In this section we consider the transitions d — > dj(g) and u —> uj(g) (denoted by d^f 1 
and du^) which are the building blocks for the electric dipole moment d n of the neutron. 
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As we are only interested in a rough estimate of d n , we do not include QCD corrections 
to these building blocks. In this approximation the latter can be described by the effective 
Hamiltonian 



K/f U = C R m d dv,vP R d + cf m d da, u P L d + 

c dd g m d da^P R T a d G a ^ + c d L d g m d da^P L T a d G a '^ + (d -> u) . (101) 

The effective Wilson coefficients cj^^ and c d ^ u g again receive neutral and charged Higgs 



contributions. The neutral contributions of the Wilson coefficients (involved in d d f ) read 

3 



cf° = E 

k,j=i 



1927r 2 m 2 „ o 



r LRH°* LRH° k LRH° k * LRH° k 
1 ddj 1 ddj + 1 djd 1 djd 



m d] LRHl LRHl 

m d dd > d ' d 



9 + 61n 



m 



(102) 



dd,H 

2 R,g 



E 



1927r 2 m^o 



r LRH°* r LRH° v LRH°± v LRH° k 



m dj LRH° k LRH° k 

m7 ^ 



9 + 61n 



(103) 



, <M,H£ _ dd,H°* 
and C L,( 9 ) - C /?,( 9 ) 

(involved in dj^) read 

3 



The charged Higgs penguin contributions to the Wilson coefficients 



',H d 



R 



E 



167T 2 m 2 

j=l u i 



r 



LRH ± *- r LRH ± riO 



da , 
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+r 



LRH ± - r LRH ±rrlu j riO 



da. 



m d 



XY 



m 



H+ 



m 



H+ 



r 
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1 u,d L/ ' 



7,YY 



m 



H+ 



(104) 



E 



167T 2 m 2 

j = l U 3 



r LRH ± * r LRH ± r^0 
1 du. 



du-j 



a 



8,YY 



i - r LRH ± - r LRH ± rrlu j r<0 
+ l dm L --* U < 



m d 



8,XY 



, dd,#± _ dd,H±* 
allClC L,( 5 ) - C R,( 9 ) • 

The analogous expressions for c^', g s 
given as 

6 



rn 



m; 



i V LRH ± * V LRH ± ^0 



rn 
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which are involved in the expressions of di 9) arc 



2 R 



E 



167r 2 m 2 
j,fc=i u j 



J- UUi J- tllli >— 's 
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J- u,n J- 11,-u t — 's 



m; 



rn 
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LRH°LRH°m Uj 
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rn 
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(106) 
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uu,H<j> 

z R,g 



E 



167r 2 m 2 
j,k=i u j 



-p LRH k*-p LRH k r^O 
J- uuj J- uuj *-^g 



YY I 2 

n k . 



J- UjU J- UjU ^8 



LRH° LRH° k m U] 
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(107) 



C R 



E T 



— e 



1527T 2 m 



^LRH ± - r LRH ± md i 



udj 



m,, 



12 In 



<ijM djM 



'1081 



iu,H% 



Again, we have c L , g s 
given in Eq. (1831) . 



E 



192tt 2 m 2 , 

3=1 H+ 



^LRH ± *- r LRH ± , r LRH ± *- r LRH ± 



udj 



+ r 



-r 



2 

LRH±T,LRH± m dj ( n , c f m rf 3 



9 + 61n 



and c r 



-i,(s) 



UU,H^ : 



(109) 



The loop functions Cysxyyy(?/i) are 



G. Leptonic decays of neutral mesons 

The effective Hamiltonian H e ff which includes the full set of operators for the general 
decays PS(q~fqi) — > (PS' refers to the pseudo-scalar meson) reads 



njAF=l 
H eff 



[C*'*0*f* + CT'OT + C^OT + ('!■'()",;'' + Primed] + h.c. 

(110) 



where the operators (together with their primed counterparts) are defined as 
= (qf^PRqi) (IbY£a) , Of* = (#7 M P*ft) (£ B jVa) 

oT = (QfPLqi) Wa) , oT = (qfPLqi) (Wa) , 

= (q f P Rqi ) (I B £ A ) , 0'¥ qi = (q f P Rqi ) (£ Bl5 £ A ) . 

Making use of the hadronic matrix elements 



(0|g/7«75<?i|PS) = ifpsPps, 

M 2 

(0\q n5qi \PS) = -if F 



I PS" 



m qf +m q .) 



Ill) 



;ii2) 
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one obtains the branching ratio 

14: ]\ #4 



B [PS(q f qJ -+ Z\r B ] 



f(x\,x%) Mpsftg {m A +m lB f 7 



X 



+ 



{m qf + m q .) {mi A + mi B ) V 

+ (mi A -m lB ) r mi _ 
(m qf + m qi )(m lA + m lB ) (m u + m lB ) \ v v 



x [1 - (x A - x B ) 2 } 



x [1 - (x A + X B ) 



where the function f(xi,Xj) and the ratio defined as 
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113) 



f{Xi,Xj) = yjl - 2(Xi + Xj) + (Xi - ~Xj)' 



.1 j 



M, 



PS 



1. Wilson coefficients 



Tree level neutral Higgs contributions to PS(qfqi) — > t\t B in the 2HDM of type III 

The non-vanishing Wilson coefficients of the operators in Eq. (IllOp induced through 
tree-level neutral Higgs (H® = H°, h°, A ) exchange read 



C p fqi 



Up 



IT 



X 1 



2G 2 F M 2 W £ m\ 



v LRH° k v RLH° k \ r RLH° 
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E( -p LRH k I -p RLH k\ -pLRH% 
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2G F M W k=l nip-o 
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i A ) 1 IfH 



(114) 



Loop-induced charged Higgs contributions to B s — > fi + fi in the 2HDM of type II 

As mentioned earlier, we also include in our analysis the 2HDM of type II loop-induced 
charged Higgs contributions to B s — > fi + fi~ from Ref. |52|: 



bs _ rn b V* b V ts m p 2 log (m 2 H /m 2 ) 
2M 2 W P m 2 H /m 2 - 1 



(115) 



where mj and m t are understood to be running masses evaluated at the matching 
scale. 
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H. Flavour-changing lepton decays 



The general expressions for the branching ratios of r — > e and r — > \i 

have the form 
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H° h 



Note that the (not explicitly denoted) sum over the Higgses must be performed before taking 
the various absolute values in Eq. (j!16p . 

I. Input parameters 

In this section we list our input parameters in tabular form. 



Parameter 


Value (GeV) 


m u (2 GeV) 


0.00219 ± 0.00015 [148] 


m d {2 GeV) 


0.00467 ± 0.00020 [148] 


m s (2 GeV) 


0.095 ± 0.006 [148] 


m c {m c ) 


1.28 ±0.04 [149] 


m b (m b ) 


4.243 ± 0.043 [88] 


m t (m t ) 


165.80 ±0.54 ±0.72 [14] 



Parameter 


Value 


M w 


80.40 GeV 


M z 


91.19 GeV 


a s (M z ) 


0.119 


G F 


1.16637 x 10~ 5 GeV" 2 


"em 


137 


V 


174.10 GeV 



TABLE IX: Left: Input values for the quark masses used in our article. In the numerical analysis, 
we used the NNLO expressions in a s for the running (see for example Ref. 1501 ]) in order to obtain 
the quark-mass values at higher scales. Right: Electroweak parameters and the strong coupling 
constant used in our analysis. Concerning the running of a s we used NNLO expressions (given for 
example in Ref. [711] ). 
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Parameter 


Value 


Ib s 1 !b 


1.221 ±0.010 ±0.033 [14] 


fo 


218.9 ± 11.3 MeV [1371 


fD a 


249 ± 2 ± 5 MeV [14] 


fD a 1 Id 


1.188 ±0.025 [1371 


fx 


156.3 ± 0.3 ± 1.9 MeV [14] 


fx/ fir 


1.193 ±0.005 [1481 



Meson massses 


Values (GeV) 


m B ±(B0) 


5.279 


™B a 


5.367 


m D ± {D o) 


1.870 (1.865) 


m D a 


1.969 


m K±(K°) 


0.494 (0.498) 


"V±(7r°) 


0.140 (0.135) 



TABLE X: Left: Values for de c ay constants of Re f. l 4| obtained by avera ging the lattice results 



of Ref. 



137], 



151 



151 1153 . 1154 1 1551 . 1 150 . Il57l . 



1581 . 1 1591 . 116CL 



161 



162J, [l63j. Right: Meson masses 



according to the particle data group (see online update of Ref. 711]). 
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